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ABSTRACT 



We have investigated the distribution and evolution of ~3 100 intergalactic neutral hydrogen (H i) absorbers with H i column densities 
log A'hi = [12.75, 17.0] at 1.9<z<3.2, using 18 high resolution, high signal-to-noise quasar spectra obtained from the ESO VLT/UVES 
archive. We used two sets of Voigt profile fitting analysis, one including all the available high-order Lyman lines to obtain reliable H i 
column densities of saturated lines, and another using only the Lyo' transition. There is no significant difference between the Lycf-only 
fit and the high-order Lyman fit results. Combining our Lya-only fit results at 1.7 < z < 3.6 with high-quality literature data, the 
mean number density at < z < 4 is not well described by a single power law and strongly suggests that its evolution slows down 
at z < 1.5 at the high and low column density ranges. We also divided our entire Hi absorbers at 1.9 < z < 3.2 into two samples, 
the unenriched forest and the Civ-enriched forest, depending on whether Hi lines are associated with Civ at logA'civ ^ 12.2 within 
a given velocity range. The entire H i column density distribution function (CDDF) can be described as the combination of these 
two well-characterised populations which overlap at log A'hi ~ 15. At log A'hi ^ 15, the unenriched forest dominates, showing a 
similar power-law distribution to the entire forest. The C iv-enriched forest dominates at log A'hi ^ 15, with its distribution function 
as oc A'jjj' "''. However, it starts to flatten out at lower A'hi, since the enriched forest fraction decreases with decreasing A'hi. The 
deviation from the power law at log A'hi = [14, 17] shown in the CDDF for the entire H i sample is a result of combining two different 
Hi populations with a different CDDF shape. The total Hi mass density relative to the critical density is £1hi ~ 1-6 x lO^^h^' , where 
the enriched forest accounts for ~- 40% of Hhi- 
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1. Introduction 

The resonant Lya absorption by neutral hydrogen (H i) in the 
warm (~ 10'* K) photoionised intergalactic medium (IGM) pro- 
duces rich absorption features blueward of the Lye emission line 
in high-redshift quasar spectra known as the Lya forest. The Lyo' 
forest contains ~ 90% of the baryonic matter at z ~ 3 and can be 
observed in a wide range of redshifts up to z~6. Gas-dynamical 
simulations and semi-analytic models have been very successful 
at explaining the observed properties of the Lyo- forest mainly 
at low Hi column densities A^hi ^ 10'^ cm"^. These models 
have shown that the Lya forest arises by mildly non-linear den- 
sity fluctuations in the low-density Hi gas, which follows the 
underlying dark matter distribution on large scales. This inter- 
pretation also predicts that the Lya forest provides powerful ob- 
servational constraints on the distribution and evolution of the 
baryonic matter in the Universe, hence the evolut ion of galax- 
ies and the large-scale structure dCen et al. |[T994l IR auch et al. 



Aguirreet all 1200 ll: ISchave et al] l2003t lOppenheimer & Pavel 
20061k 



1997HTh eunseta l.|[T99 
Schavd I2OOI; Ki m et al 



Dave et al.lll999tlSchave et al.i,2000h: 
20021) . In addition, the discovery of 



triply ionised carbon (C iv) associated with some of the forest ab- 
sorbers suggests that the forest metal abundances can be utilised 
to probe early generations of star formation and the feedback 
between high-redshift galaxies and the surrounding IGM from 
which galaxies formed (Cowie et al. 1995; Dave et al.. .1998: 
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* The data used in this study are taken from the ESO archive for the 
UVES at the VLT, ESO, Paranal, Chile. 



The physics of the Lya forest is mainly governed by three 
competing processes, the Hubble expansion, the gravitational 
growth and the ionizing UV background radiation. The Hubble 
expansion which causes the gas to cool adiabatically and the 
gravitational growth are fairly well-constrained by the cosmo- 
logical parameters and the primordial power spectrum from the 
latest WMAP observations (IJarosik et all 120 11 ). On the other 
hand, the ionizing UV background radiation controls the pho- 
toionisation heating and the gas ionisation fraction, thus deter- 
mining the fraction of the observable H i gas compared to the 
unobservable Hii gas. The UV background is assumed to be 
provided primarily by quasars and in some degree also by star- 
forming gala xies ( Shaplev et al.l l2006t ISiana et aTl 12010 ) and 
Lyof emitters dlwata et aI.ll2009l) ~However. the intensity/spectral 
shape of the UV background and the relative contribution from 
quasars and g alaxies as a function of redshift a re not well con- 
sti-ained (Bol ton et al.ll2005t Ipaucher-Giguere et al. 2008). One 
of the common methods to measure the UV background and it s 
evolution is the quasar proximity effect dDaU'Aglio et alj|2008l) . 
Unfortunately, measurements of the UV background through the 
proximity effect are biased by the large scale density distribution 
around t he quasars which cann ot be easily quantified observa- 
tionally (Parti et aI.ll20Toll201 Ih . 

Two commonly explored quantities to constrain the proper- 
ties of the Lya forest are the number of absorbers for a given 
H I column density range per unit redshift, d«/dz, and the differ- 
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ential column density distribution function (CDDF, the number 
of absorbers per unit absorption path length and per unit col- 
umn density, an analogue to the galaxy luminosity function). 
Compared with simulations, detailed structures seen in an over- 
all power-law-like CDDF (oc A^^j) such as a flattening or a steep- 
ening at different column density ranges const rain various for- 
est physical a nd galactic feedback processes (Altay et al."2011; 
iDave etalJ26 lO). The CDDF is also one of the main observables 
required in calculating the mass density relative to the critical 
density contributed by the forest (ISchavell200lh . The shape of 



the CDDF at lower A'hi < 10 



12.5-12.7 



cm s (a typical detection 



limit for most available high-quality data) is of particular impor- 
tance, since the lower A^hi absorbers are much more numerous 
than higher A^hi absorbers, thus they can trace a significant frac- 
tion of baryons, depending on the steepness of the CCDF at the 
low A^H I limit. 

On the other hand, dn/dz provides an additional way to study 
the UV background radiation and its evolution. The gas den- 
sity decreases with decreasing redshift due to the Hubble ex- 
pansion. A lower gas density results in a strong reduction of the 
recombination rate, allowing the gas to settle in to a photoion- 
isation equilibrium with a higher ionisation fraction. With the 
non-decreasing background radiation, this causes a steep number 
density evolution. However, the decrease of the quasar number 
density at z < 2.5 also decreases the available ionising photons 
(ISilverman et al.l 12005). This changes the ionisation fraction in 
the gas and also counteracts the gas density decrease, and hence 
slows down the number density evo lution (iTheuns et al. I [19981 
iDave et alJ[T999HBianchi etani200lh . 

The result from the HST/FOS Quasar Absorption Line Key 
project s hows such a slow cha nge in the d«/dz evolution at 
z < 1.5 ("Wevma nn et al.l [T998h , compared to a much steepe r 
d«/dz evolution shown at z > 2 (iKim et alJll9"97l l200lL l2002h 
Cosmic variance also seems to increase at lower z dKim et alj 
l2002h . Unfortunately, recent work based on better-quality HST 
data at z < 1.5 (or the observed Hi Lya at < 3050 A) have 
shown rather ambig uous dw/dz results with a large scatter along 
different sight Unes (iJanknecht et al.l l2006t iLehner et all l2007t 
IWiUiger et al.ll2()10h . The only certain observational fact is that 
all of these newer studies show a factor of ~ 2 - 3 lower number 
densities than the Weymann et al. values at z< 1.5. Considering 
a lack of results from good-quality data at 1 <z < 1 .5 in the litera- 
ture, the redshift evolution of dn/dz can be considered as a single 
power law without any abrupt change in dn/dz at 0<z< 3.5. 

Here we present an in-depth Voigt profile fitting analysis of 
18 high resolution (R ~ 45 000), high signal-to-noise (~ 35-50 
per pixel) quasar spectra obtained with the UVES (Ultra-violet 
Visible Echelle Spectrograph) on the VUT, covering the Lya for- 
est at 1.9 < Zforest < 3.2. Our main scientific aims are to derive 
the redshift evolution of the absorber number density and the 
column density distribution function from a large and homoge- 
neous set of data available at z > 2, since most previous high- 
quality forest studies at z > 2 have been based on less than 5 
sightlines. Even with few sightlines, the statistics for the weak 
forest lines is robust due to the large number of weak absorbers 
with A^Hi = lO'^"'^ cm"^ (about 150 absorbers at z~2.5 per sight- 
line, i.e. in the wavelength range between the quasar's Lya and 
LyyS emission lines). However, for the stronger forest systems 
with A^Hi ^ 10'^ cm"^, more sightlines are required since there 
are only about 10 absorbers per sightline at z~2.5. Cosmic vari- 
ance also plays an imp ortant role at low er redshifts, especially 
for stronger absorbers (iKim et al.ll2002h . Therefore, increasing 



the sample size at z ~ 2 is critical in addressing the A^h i evolution 
for the Lya forest. 

In addition to the increased sample size, we have improved 
previous results in two ways. First, most previous studies on the 
forest from ground-based observations at z > 1.7 have been 
based on the LyfT-only profile fitting analysis. This approach 
does not provide a reliable A^hi for saturated lines, A^hi ^ lO'''^ 
cm"^ for the present UVES data. To derive a more reliable A^hi 
of saturated lines, we have included all the available high-order 
Lyman series in this study. 

Second, strong evidence have been accumulated in recent 
studies that metals associated with the high-redshift Lyo' for- 
est are within ~ 100 kpc of galaxies as in the circum-galactic 
medium rather than in the intergalactic space far away from 
galaxies dAdelberger et al.ll2005l: ISteidel et al.ll2010l:lRudie et"al] 
201:^. This implies that the Hi absorbers containing metals 
might show different properties than the ones without detectable 
metals. Taking Civ as a metal proxy, we have divided our data 
into two samples, one with C iv (the C iv-enriched forest) and 
another without C iv (the unenriched forest), in order to test this 
scenario of the circum-galactic medium. Since our study lacks 
the imaging survey around the quasar targets, we cannot claim 
that the Civ-enriched forest is indeed located within ~ 100 kpc 
from a nearby galaxy. However, this study provides complemen- 
tary results to galaxy-absorber connection studies at high red- 
shifts (Steidel et al. 2010; Rudie et al. 2012). 

This study is also very timely since the Cosmic Origins 
Spectrograph (COS), a high-sensitivity FUV spectrograph on- 
board HST has s tarted to produce many high -quality quasar 
spectra at z < 1 dGreen et al.ll2012l: ISavage et al.l l2012). These 
COS quasar observations have opened a new tool to study the 
low-z Lya forest. Combined with results at high redshifts such 
as our study, COS observations will make it possible to charac- 
terise the d«/dz evolution at < z < 3.5 in a more robust way, 
thus a stringent constraint on the UV background evolution. 

This paper is organised as follows. Section |2] describes 
the analysed data and two Voigt profile fitting methods. 
Comparisons with previous studies based on the Lyor-only fit 
are shown in Section |3] The analysis based on the high-order 
Lyman fit is presented in Section|4] Column density distribution 
and evolution of the Lya forest containing C iv are presented in 
Section|5] Finally, we discuss and summarise the main results in 
Section|6] All the results on the number density and the differen- 
tial column density distribution from our analysis are tabulated 
in Appendix A. Throughout this study, the cosmological param- 
eters are assumed to be the matter density 0,„ - 0.3, the cos- 
mological constant Qa = 0.7 and the current Hubble constant 
H() - lOOh km s ' Mpc"' with h - 0.7, which is in concordance 
with latest WMAP measurements (Jarosik et al. 201 1). The log- 
arithm A^hi is defined as logA^ei - log(A^Hi/l cm"^). 

2. Data and Voigt profile fitting 

Table [T] lists the properties of the 18 high-redshift quasars anal- 
ysed in this study. The redshift quoted in Column 2 is mea- 
sured from the observed Lya emission line of the quasars. 
Note that the redshift based on the emission lines is known 
to be under-estimated compared to the one measured from 
other quasar emission li nes such as Civ dTvtler&Fanlll99a 
IVanden Berk et al.l 1200 lb . The spect r um of Ql 101 -264 is the 
same one as analysed in Kim e t aP d2002h . while the rest of 
spectra are from Kimetal. (2007|). The raw spectra were ob- 
tained from the ESO VLT/UVES archive and were reduced with 
the UVES pipeline. All of these spectra have a resolution of 
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Table 1. Analysed quasars 



Quasar 




~Lya 




"Ly.high-order 


Excluded 


S/N per pixel^" 


LL (A) 


notes 


Q0055-269 


3.655' 


2.936-3.605 




2.936-3.205 


not used 


50-80 [..., ...] 


2288 




PKS2126-158 


3.279 


2.815-3.205 




2.815-3.205 




50-200 [..., 170] 


3457 


two sub-ULAs at z =2.768 & 2.638 


Q0420-388 


3.116' 


2.480-3.038 


2.665-3.038 


2.665-3.038 


2.607-2.670 


100-140 [..., 120] 


3754 


a sub-DLA at z =3.087 


HE0940-1050 


3.078 


2.452-3.006 






2.714-2.778 


50-130 [..., 115] 


< 3200 




HE2347^342 


2.874' 


2.336-2.819 






2.708-2.773 


100-160 [105, 120] 


< 1160 


multiple associated systems 


Q0002-422 


2.767 


2.209-2.705 








60-70 [140. 170] 


3025 




PKS0329-255 


2.704' 


2.138-2.651 








40-60 [90. 90] 


3157 


an associated system at 4513.7 A 


Q0453^23 


2.658' 


2.359-2.588 
2.091-2.217 








90-100 [130, 120] 


3022 


a sub-DLA at z =2.305 


HE 1347-2457 


2.609' 


2.048-2.553 








85-100 [115, 120] 


2237 




Q0329-385 


2.434 


1.902-2.377 








50-55 [(50, 85), ...] 


< 3050 




HE22 17-28 18 


2.413 


1.886-2.365 


1.971-2.365 


1.970-2.365 




65-120 [125, ...] 


2471 




QO 109-35 18 


2.405 


1.905-2.348 


1.974-2.348 


1.968-2.348 




60-80 [(80, 140), ...] 


2163 




HEl 122-1648 


2.404 


1.891-2.358 


1.974-2.358 


1.900-2.358 




70-170 [240, ...] 


< 1629 




J2233-606 


2.250 


1.756-2.197 


1.980-2.201 


1.970-2.201 


not used 


30-50 [...,...] 


< 1750 




PKS0237-23 


2.223' 


1.765-2.179 


1.974-2.179 


1.961-2.179 




75-110 [(130, 200), ...] 


< 3050 


a sub-DLA at 2=1.673 


PKS1448-232S 


2.219 


1.719-2.175 


1.974-2.175 


1.953-2.175 




30-90 [(70, 120), ...] 


< 3050 




QO 122-380 


2.193 


1.700-2.141 


1.977-2.141 


1.977-2.141 




30-80 [(85, 130), ...] 


< 3052 




QllOl-264 


2.141 


1.882-2.097 


1.970-2.097 


1.900-2.097 




80-110 [90, ...] 


2597 


a sub-DLA at z =1.839 



" The redshift range of the Lya forest region analysed for the number density evolution in the Lya-only fit. For the differential column density 
evolution, we used the redshift range listed in Column 5. 

The redshift range for which the high-order Lyman fit can be performed is listed only when it is different from the Lya-only fit region. 
The redshift range of the Lyo- forest region analysed for the high-order Lyman fit is listed only when it is different from the Lycf-only fit region. 
The redshift range excluded for the C iv-enriched H i forest analysis in Section 5 due to the lack of the coverage. Q0055-269 and J2233-606 
are excluded due to their lower S/N in the C iv region. No entries mean that the analyzed zciv is the same as the one in Column 5. 
° The number outside the bracket is a S/N of the H i forest region. The first number inside the bracket is a typical S/N of the C iv region at 
1.9 < z < 2.4, while the second is for 2.4 < z < 3.2. The dotted entries indicate that no C iv forest region is available for a given redshift range. 
Two numbers inside the parentheses indicate the S/N of the Civ region at 1.9 < z < 2. 1 and 2.1 < z < 2.4, respectively, due to the dichroic 
setting toward some sightlines. 

' Due to the absoiption systems at the peak of the Lyo- emission line or to the non- single-peak emission line, the measurement is uncertain. 
8 Part of the continuum uncertainties toward shorter wavelengths is due to the local, non-smooth continuum shape, partly due to a lower S/N 
(-30-35). 



R~45 000 and heliocentric, vacuum-corrected wavelengths. The 
spectrum is sampled at 0.05A. A typical signal-to-noise ratio 
(S/N) in the Lya forest region is 35-50 per pixel (here after all the 
S/N ra tios a re given as per pix el). Readers can refer to lKim et alJ 
(12004') and ' Kim et alj ( 12007 ^ for the details of the data reduc- 
tion. To avoid the proximity effect, the region of 4,000 km s ' 
blueward of the quasar's Lya emission was excluded. 

In order to obtain the absorption line parameters (the red- 
shift z, the column density N in cm"^ and the Doppler parameter 
b in km s"'), we have performed a Voigt profile fitting analysis 
using VPFIlfl Details can be found in the documentation pro- 
vided with the s oftware, Carswell, Schaye & Kim (2002) and 
iKim et al. I (l200l . Here, we only give a brief description of the 
fitting procedure. 

First, a localised initial continuum of each spectrum was 
defined using the CONTINUUM/ECHELLE command in IRAF. 
Second, we searched for metal lines in the entire spectrum, 
starting from the longest wavelengths toward the shorter wave- 
lengths. We first fitted all the identified metal lines. When metal 
lines were embedded in the H i forest regions, the H i absorption 
lines blended with metals are also included in the fit. Sometimes 
the simultaneous fitting of different transitions by the same ion 
reveals that the initial continuum needs to be adjusted to obtain 
acceptable ion ratios. In this case, we adjusted the initial con- 
tinuum accordingly. The rest of the absorption features were as- 
sumed to be H I. 

After fitting metal lines, we have fitted the entire spectrum 
including all the available higher-order Lyman series in the 
UVES spectra. This is absolutely necessary to obtain reliable 

' Carswell et al.: http://www.ast.cam.ac.uk/~rfc/vpfit.htmll 



Nui for saturated lines, as our study deals with saturated lines 
and relies on line counting. A typical z ~ 3 IGM absoiption 
feature having b ~ 30 km s"' starts to saturate around Nhi > 
10''*^cm"2 at the UVES resolution and S/N. Unfortunately, Nm 
and b values of saturated lines are not well constrained. In order 
to derive reliable Nh i and b values, higher-order Lyman series, 
such as Ly^S and Lyy, have to be included in the fit, as higher- 
order Lyman series have smaller oscillator strengths and start to 
saturate at much larger Nhi- 

During this process, another small amount of continuum re- 
adjustment was often required to achieve a satisfactory fit, i.e. a 
reduced value of ~ 1.2. With this re-adjusted continuum, we 
re-fitted the entire spectrum. This iteration process of continuum 
re-adjustments and re-fitting was then repeated several times un- 
til satisfactory fitting parameters were obtained. This produces 
the final set of fitted parameters for each component of the high- 
order Lyman fit analysis. 

In addition to this high-order fit, we have also performed the 
same analysis using only the Lya transition region, i.e. the wave- 
length range above the rest-frame Ly/? and below the proximity 
effect zone. This additional fitting analysis was done, since most 
previous studies on the IGM Nhi analysis based on Voigt pro- 
file fitting utilised only the Lya region. For the Lya-only fit, we 
kept the same continuum used in the high-order Lyman fitting 
process. In principle, the difference between two sets of fitted 
parameters occurs only in the regions where saturated absorp- 
tion features are included. In both fitting analyses, we did not tie 
the fitting parameters for different ions. 

Fig. [T] shows the numbers of absorption lines as a func- 
tion of Nhi for both fitting analyses at the two redshift ranges, 
2.2 < z < 2.6 and 2.8 < z < 3.2, in order to illustrate the differ- 
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Fig. 1. Numbers of absorption lines as a function of A^hi at 2.2 <z<2.6 and 2.8 <z<3.2. The Lyc-only fits are shown as soHd lines, 
while the high-order Lyman fits are marked as dashed lines. Solid errors indicate the lo" Poisson errors of the Lyc-only fits. 



ences at high and low redshifts. The differences between the two 
samples occurs mostly at A^hi ^ lO''*'^ cm"^ and at A^hi ^ 10'^ 
cm"^. This difference in the line numbers at A^hi ^ 10'^'^ cm"^ 
seems to be stronger at 2.8 < z < 3.2, although it is still within 
2cr Poisson errors. The line numbers at A^hi ^ 10'^ cm"^ are 
more susceptible to the incompleteness which depends on the 
local S/N than the difference between the two fitting methods. 
The difference at other column density ranges is smaller, which 
in turn leads us to expect that there is no significant difference 
between the Lya-only fit and the high-order Lyman fit. 

We restrict our present analysis to log A^hi = [12.75, 17] at all 
redshifts. As clearly seen in Fig.[T] the incompletenes s becomes 
quite severe for logA^Hi < 12.5 and redshifts z > 3 dKim et alj 
119971 l2002h . Therefore, the lower A^hi limit was chosen to be 
logA^Hi = 12.75. We chose logA^Hi = 17 as the upper A^hi limit 
since we wanted to analyze only the Lya forest whose traditional 
definition is an absorber with log A^hi < 17.2 (above which it is 
referred to as a Lyman limit system (Tytler 1982)). Additionally, 
absorbers at log A^h i > 17 are very rare (Fig. [1]. 

Note that the availability of the high-order Lyman series de- 
pends on the redshift of the quasar and whether the sightline con- 
tains a Lyman limit system. In addition, the amount of blending 
affects whether a reliable column density can be measured. At 
high redshifts Zem > 3, line blending becomes severe. However, 
most UVES spectra also covers down to 3050 A where Lyman 
lines higher than Ly6 are available. On the other hand, at Zem < 
2.5 the available high-order Lyman lines are rather limited, with 
mostly Lyy6 and Lyy available. However, line blending is less 
problematic than at higher redshifts. We have generated tens of 
saturated artificial absorption lines and fitted them including and 
excluding high-order Lyman lines. These simulations show that 
unblended absorption features at A^hi ^ 10'^ cm"^ can be rea- 
sonably well constrained with Lya and LyyS only. This indicates 
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Fig. 2. Total absorption distance X{z) covered with our sample 
of 18 high-redshift quasars. The solid line is for the Lya-only fit, 
while the dashed one is for the high-order fit. 



that our A^hi can also be considered reliable even at z < 2.5 with 
Lya and LyyS only. 

The absorption distance is obtained by integrating the 
Friedmann equation for a 0,„ - 0.3 and Qa = 0.7 universe, 
and is given by 

= J ^^{l+zfdz (1) 

(lBahcall&Peebleslll969[) . where Hq is the Hubble constant at 
z - 0. The total absorption distance X{z) covered by the spectra 
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for both Lyo'-only and high-order fits is shown in Fig.|2] The red- 
shift coverage of our sample steadily increases with decreasing 
redshift until it reaches its maximum at z~2.1. For redshifts be- 
low z < 1 .9 the coverage decreases rapidly and our sample ends 
at z = 1.7. Note that the lowest redshift possible for the high- 
order Lyman line analysis is z ~ 1.97, while the Lyo'-only fit 
analysis is possible down to z ~ 1.7. Due to the reduced redshift 
coverage in the high-order Lyman range of individual sight lines 
caused by intervening Lyman limit systems, the sample coverage 
of the high-order fit analyses is reduced between 2.4 < z < 2.7. 
At the high redshift end z > 3.22, the number of available for- 
est lines decreases and the sample consists of only one line of 
sight. The low redshift limit for the high-order fit was set to be 
the lowest redshift without any saturated lines when no Ly/? is 
available for each quasar. This criterion restricts our high-order 
Lyman fit analysis to 1.9<z<3.2. Since the redshift coverage of 
low-z quasars for the high-order fit is shorter than the one for the 
Lyg-only fit and the high-A^Hi forest clusters stronger at lower z 
(iKim et alJl2002h . the quasar-by-quasar dn/dz at z ~ 2 from the 
high-order fit analysis is expected to suffer from the low number 
statistics. 

In Table [1] Columns 3-6 summarise the redshift range used 
for the different analysis. Column 3 lists the redshift range of 
the Lya forest region analysed for the number density evolu- 
tion in the Lycc-only fit. For the differential column density evo- 
lution of the Lya-only fit, we used the redshift range listed in 
Column 5. Column 4 and Column 5 list the redshift range for 
which the high-order Lyman fit can be performed and the one 
for which the high-order Lyman fit was done, respectively. The 
region is listed only when it is different from the Lya-only fit 
region in Column 3. Since there are no strongly saturated Lya 
lines at 1.90 < z < 1.98 for some low-z quasar sightlines, we 
used a lower redshift range than the one listed in Column 4 for 
the high-order Lyman fit analysis for these sightlines. Column 6 
shows the redshift range excluded for the C iv-enriched H i study 
in Section 5. Due to the wavelength gaps caused by the UVES 
dichroic setup, the covered C iv redshift ranges are smaller than 
the Lya forest coverage listed in Column 3. The region of +200 
km s"' from the gap was excluded, and only the redshift range 
covering both C iv doublets was included in the analysis. The 
blank entries mean that the analyzed zc iv is the same as the for- 
est ZLy.high-order- Q0055-269 and J2233-606 are excluded in the 
Lya-C IV forest study due to their lower S/N in the C iv region. 

In the HE2347-4342 Lya forest region, there are very strong 
O VI absorptions mixed with the two saturated Lya absorption 
systems at 4012-4052 A dFechner et al.ll2004 . Since the fitted 
line parameters for these Lya systems cannot be well constrained 
(their corresponding hyfi is below the partial Lyman limit pro- 
duced by the z ~ 2.738 systems), we excluded this forest region 
toward HE2347-4342. In the J2233-606 sightline, there are two 
partial Lyman limit systems at 3489 A (z ~ 1.870) and 3558 A 
(z ~ 1 .926) and several high column density forest absorbers at 
3400-3650A. To derive a robust A^hi, we included the HSTjSTlS 
echelle spectrum of J2233-6043 at 2280-3 150 A. The resolution 
in this wavelength region is ~ 10 km s"' and its S/N is ~ 8 per 
pixel (0.05 A). 

Table[T]also lists the S/N of each quasar spectrum in Column 
7. The number outside the bracket is a S/N of the Hi forest 
region. The first number inside the bracket is a typical S/N 
of the Civ region at 1.9 < z < 2.4, while the second is for 
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^ The STIS spectrum is taken 

|http : //www . stsci . edu/ f tp/observing/hdf /hdf south/hdf s . html| 

dSavaglio et al.lll999h . 



from 



4.0 



Fig. 3. Number of Hi absorbers with logA^Hi = [12.75,17] 
as function of redshift in our our sample of 18 high-redshift 
quasars. The solid line is for the Lya-only fit. The dashed line 
is for the high-order fit, while the red heavy dot-dashed line is 
for the Lya-only fit for the redshift range used for the high-order 
fit. 



2.4 < z < 3.2. The dotted entries inside the bracket indicate that 
no C IV forest region is available for a given redshift range. The 
low redshift bin of the C iv forest covers the wavelength region 
where the different CCDs from two dichroic settings were used 
at ~ 4780 A(or z ~ 2.1). This leads to a much lower S/N at 
< 4780 A (z < 2.1). When the lower S/N region is larger than 
20% of the whole C iv forest range, two numbers were listed in- 
side the parentheses. The first number corresponds to the lower 
S/N at 1.9 < z < 2.1, while the second number is for the higher 
S/N at 2.1 <z<2.4. 

In addition, regions within +50 A to the center of a sub- 
damped Lya system (A^hi ^ lO'^ cm"^) are excluded, since 
they are associated directly with intervening high-z galactic 
disks/halos and could have a possible influence on the appar- 
ent line densities in the forest. The sightline toward Q0453-423 
includes a sub-DLA, which introduces a gap in the Lya red- 
shift range. All the calculations toward Q0453-423 account for 
this redshift gap correctly. However, they are plotted as a sin- 
gle data point and their plotted redshift range is the whole Lya 
redshift range without showing a gap. The sightlines toward 
PKS2126-158 and Q0420-388 also contain an intervening sub- 
DLA, which shortens the continuously available redshift cover- 
age for the high-order fit. Column 8 of Table [Tjlists the observed 
wavelength of a Lyman limit (LL, 912 A in the rest-frame wave- 
length) of each quasar, which is defined as the wavelength be- 
low whi ch the observe d flux becomes 0. The values are taken 
from Kim et al. I (I2OOI . When a Lyman limit is not detected 
within available data, it is denoted to be less than the lowest 
available wavelength. Column 9 of Table [T]notes information on 
sub-DLAs along the sightline. When a sub-damped Lya system 
exists along the sightline, we discarded 50 A centred at the sub- 
DLA each side to exclude its influence on the forest, such as a 
higher frequency or lack of higher-column density forest. The to- 
tal number of Hi lines for logA^Hi = [12.75, 17] at 1.9 < z < 3.2 
is 3077 for the high-order Lyman fit sample. The Lya-only fit 
sample has 3778 Hi lines at the total redshift range listed in the 
3rd column of Table [1] 
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In Fig. [3] the number of Hi absorbers with logA^Hi - 
[12.75, 17] from both fitting methods is shown as a function 
of redshift. The number of absorbers obtained from each fit- 
ting analysis is roughly proportional to the absorption distance 
coverage. Therefore, our sample shows the highest H i absorber 
numbers around redshift z ~ 2 for each fitting analysis, where 
the sample absorption distance coverage also reaches its maxi- 
mum. Sometimes the high-order fit analysis (dashed line) reveals 
a slightly higher number of absorbers between 2 < z < 3. This is 
because what appear to be single saturated Lya lines may have 
more than one component present in the corresponding higher 
order Lyman lines. At z ~ 2, the number of the Lya-only-fit ab- 
sorbers (heavy dot-dashed line) is slightly larger than the high- 
order-fit absorbers. This is caused by the fact that some simple 
saturated lines with log A^hi < 17 in the Lya-only fit analysis are 
actually absorbers with log A^hi > 17 in the high-order fit analy- 
sis. Since the Lya-only fit gives a lower A^h i limit for a saturated 
line, these lines are included in the Lya-only fit sample, but ex- 
cluded in the high-order fit sample in Fig.|3] 

3. Comparison with previous studies using Lya 
only 

In Sect. |2] we have shown that including higher order transitions 
in the fitting process slightly alters the column density statistics 
at logA^Hi > 15.0. In order to compare our quasar sample with 
previous studies based only on the Lya transition, we briefly 
present the column density distribution and evolution derived 
from the Lyc-only fit in this section. A large redshift coverage 
is very important in the study of the absorber number density. 
Therefore we used all Lya lines found in the whole available 
Lya redshift ranges listed in Column 3 of Table[T]in this section. 
On the other hand, the diff'erential density distribution function 
is not sensitive to a large redshift coverage. Thus, only the Lya 
lines at 1.9 < z < 3.2 are analysed for the distribution function 
study. A detailed analysis using the high-order fit is presented 
in Section |4] All the results from this section are tabulated in 
Appendix A. 

3.1. Absorber number density evolution dn/dz 

The absorber number density n(z) is measured by counting the 
number of H i absorption lines for a given column density range 
for each line of sight. The line count n is then divided by the 
covered redshift range Az to obtain dn/dz. If forest absorbers 
have a constant size and a constant comoving number density, 
its number density evolution due to the Hubble expansion can 
be described as 



^ = nR^NocH(zr\l + z)\ 
dz 



(2) 



where R is the size of an absorber, A^o is t he local comoving 
number density and c is the speed of Ught (iBahcall & Peebles 
[19691) . For our assumed cosmology, Eq.|2]becomes 



dn 



(3) 



dz V0-3(l +z)3 -hO.7 

At 1 < z < 4.5, Eq.[3]has an asymptotic behaviour of dn/dz oc 
(1 + z)~°'', while at z < 1 it becomes dn/dz « (1 -h z)~' '^. 
For higher redshifts the asymptotic behaviour becomes dn/dz oc 
(1 -H z)'^^. Any differences in the observed exponent from what 
is expected from Eq. |3]indicate that the absorber size or/and the 
comoving density are not constant. 
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Fig. 4. The number density evolution of the Lya forest in the col- 
umn density range log A^h i = [14, 17] of the Lya-only fits. Black 
filled circles show results from our data set, which is tabulated 
in Table lA.fl Other data points indicate various results obtained 
from the literature. The vertical error bars give the 1 cr Poisson 
error, while the x-axis error bars show the redshift range covered 
by each sightline. The solid line shows the fit to our data only. 
Dashed line is the result including the literature data for z > 1 

S log(l -Hz) > 0.3). The dotted line gives the fit given in lKim et alJ 
2002h . The green dot-dashed curve shows the predicted dn/dz 
evolut ion based on a quasar-only UV background by Dave et aT] 
(Il999h . The red dotted and the blue dot-dot-dot-dashed curves at 
z < 2 illustrate the predicted dn/dz based on momentum-driven 
wind and no-wind mo dels with a UV ba ckground by quasars and 
galaxies, respectively (iDave et al.ll2010l) . 



Empirically, dn/dz is described as dn/dz = A(l -i- zY. It has 
been known that dn/dz evolves more rapidly at higher column 
densities. At z > 1.5, a y ~ 2.9 is found fo r A^hi = lO'"^"'^ 

(iKiin etal 



cm 



-2 



and y 



z > a y 
1.4 for A^H 



10' 



cm 



an i2002h . 

At z < 1 .5, IWevmann etall (1 19981) found y ~ 0. 16 and A ~ 35 for 
absorbers with a rest-frame equivalent width greater than 0.24 A 
from HST/FOS data. Later studies on dn/dz based on the pro- 
file fitting or curve of growth analysis using better-quality data 
from HST/STIS and HST/G HRS show a factor of -2-3 lower 
dn/dz than the one found bv Wevmann et al.] (Il998h . These stud- 
ies also show a l arger scatt er in dn/dz at z < 0.2 with A ~ 5- 
22 (iLehner et al.l l2007: Williger et al.|l2010). Part of this scatter 
is thought to be caused by inhomogeneous data quality, anal- 
ysis methods, and cosmic variance. Unfortunately high-quality 
data lack a complete z coverage at z < 1.5, missing mostly at 
0.4 < z < 1.0. Keep in mind that the EOS result and most avail- 
able ground-based results at z > 1 .5 in the literature are based 
on the Lya Unes only, while most space-based results at z < 1 .5 
are using the available high-order Lyman series. Therefore, it is 
not possible to derive a robust power-law slope y of dn/dz at 
< z < 3.5. Strictly speaking, a fair comparison should be made 
on the data with similar qualities and uniform analyses. 

The number density evolution is illustrated in Figs. |4]|5] and 
|6]for two different column density ranges: logA^Hi = [14, 17], 
and [13.1,14]. Data compiled from the literature are indi- 
cated in t he fig ures: Hu et al. ( 1995), Lu et al. ( 1996), Kim et slU 
(Il997h . liciAman & Tvtler. (.19971 .Wevmann et aL (,199^ 
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Savaglio et al l (Il999h iKim ^ ll2001h. ISembach et al. I (12004 
Williger et al.' ('2006V 'Aracil et al.' ('20061^1 'Janknec ht et al.1 
(2006fl Lehneretal. (2007) and Williger et al. (2010). To be 
consistent with our definition of the proximity effect zone, we 
applied the same 4 000 km s"' exclusion within the quasar's Lya 
emission line for all the literature data, whenever the line lists 
from the literature include all the Lya lines below the Lyo- emis- 
sion line of the quasar. When the published line lists are only 
for the shorter wavelength region than the entire, available for- 
est regi on outside the 40 00 km s ' proximity zone, such as the 
ones bv iHu et al.1 ([1991 . no such an exclusion is required. We 
used all the reported H i lines in the literature mentioned above, 
without any pre-selection imposed on A^h t or b parameters. The 
latest study on the low-redshift IGM bv I Williger et all (1201 Ol) 
found that the number density from the HST/ STIS results is a 
factor of 2-3 lower than the HST/FOS results bv lWevmann et al.) 
(119981). They apphed the same selection criteria on H i absorbers 
used bv iLehner et al.l (l2007h . i.e. measurement errors less than 
40% and b <4Q km s"' . As H i absorbers tend t o have a larger b 
parameter at lower redshift (iLehner et al.ll200'7l ) and larger mea- 
surement errors in general, selecting H i absorbers at > 40 km 
s"' has a larger impact on dn / dz at lower redshift. In addition, as 
the HSTfFOS results are based on the H i sample without any im- 
posed selection criteria, using the full H i lines provides a more 
straightforward comparison to the HSTfFOS result. 

We have performed a linear regression to our data in log- 
arithmic space for the various column de nsity bins, using the 
maxim um likelihood method described in iRiplev & ThompsonI 
(Il987h . This method accounts for the uncertainties in the num- 
ber density and incorporates the weighting using the uncertain- 
ties. Errors of the fit parameters were obtained using the maxi- 
mum likelihood method. Linear regressions were once obtained 
from our data including the literature data and once without 
them. Since for redshifts z < 1 (or log(l -H z) < 0.3) the 
nu mber density evolution could remain constant with redshift, 
cf ■ IWevmann et al] (Il998h . only the literature data with redshift 
z > 1 was used for the fit. The resulting parameters are given in 
Tables 12 

Fig. ID shows the d«/dz evolution for the column density in- 
terval of log A^h i = [14, 17]. Ourresults (filled circles) agree well 
with previous findings at z > 1.5 (log(l -i- z) > 0.4), confirming 
that there is a real sightline variation in dw/dz. lKime t al. (2002) 
notes that the scatter between different sightli nes increases as 
z decr eases down to z ~ 2. In fact, the data of iJanknecht et alJ 
( 120061) at redshifts below z ~ 2 (log( 1 H-z) ~ 0.45) indicate that the 
scatter might well increase at lower z, although the errors are still 
very large to draw any firm conclusions. Considering that the 
FOS result is based on the equivalent width measurement, and 
the conversion from the equivalent width to the column density 
requires the b parameters of individual absorbers, which are ill- 
constrained at the FOS resolution, the full HST/STIS H i sample 
toward some sightlines is in good agreement with the HSTfFOS 
result (blue open triangles), although there still is a large sight- 
line variation. The full H i sample at z < 0.4 strongly supports 
the previous conclusion obtained by the HSTfFOS result, that 
dn/dz flattens out at z < 1.5. 



^ The revised hne list was used. 
The fitted hne parameters bv IJanknecht et al .1(1200 6') show many H 



lines with b < 20 km s"', about 25% of all lines. They attribute this to 
their low signal-to-noise data of less than 10 per resolution element. 
Although there are 9 sightlines analysed, one sightline has a long wave- 
length coverage from VLT/UVES and HSTfSTlS. This sightline was 
split into two data points in Figs.|4]and[5] 



Table 2. Linear regression results for d«/dz 

a) Using the Lycf-only tits and fits including literature data for z > 1 
UVES quasar by quasar quasar by quasar with lit. 



A log 


A'h, 


log A 


r 


log A 


r 


13.64- 


- 17.0 


0.78+0.14 


2.55 + 0.27 


1.13 + 0.06 


1.90 + 0.11 


14.0- 


17.0 


0.03 ± 0.20 


3.40 + 0.36 


0.72 + 0.08 


2.16 + 0.14 


14.5- 


17.0 


-0.14 ±0.26 


3.02 + 0.48 


0.46 + 0.11 


1.94 + 0.20 


13.1 - 


14.0 


146 + 0.11 


1.67+0.21 


1.52 + 0.05 


1.51+0.09 


12.75 - 


- 14.0 


1.98+0.08 


1.13 + 0.16 








b) The 


mean dn/dz using the Lyo- 


only fits at < z < 4 


A log 


A'h, 


log A 


r 






13.1 - 


14.0 


1.88 + 0.03 


0.89 + 0.06 






14.0- 


17.0 


1.03+0.07 


1.61+0.12 







c) The UVES high-order fits at 1.9 < z < 3.2 



AlogA^Hi 



quasar by quasar 
log A 7 



Mean sample 
log A 7 



12.75 - 14.0 1.85 +0.12 1.38 + 0.22 1.89 + 0.13 1.28 + 0.24 
14.0- 17.0 -0.76 + 0.29 4.91 +0.53 -0.65 +0.36 4.65 +0.66 



The linear regression to our results only (the solid line) with 
y = 3.40 + 0.36 is different at 3cr from the fit to all the available 
data at z > 1 (log(l +z) > 0.3) which yields y = 2.16 + 0.14 (the 
dashed line). This discrepancy is mainly due to the sparse data of 
our sample at higher redshift z > 3.5 (log(l + z) > 0.65) and the 
missing constraints at z < 2.0. The discrepancy is also in part 
caused by how the power-law fit is performed. Our maximum 
likelihood fit does the weighted fit. This gives a higher weight 
on higher-z data points where the Icr Poisson error is usually 
smaller The non-weighted fit for our UVES data only results in a 
steeper power-law slope, (-0.84+0.24)x(l H-z)^ °^="'-^''. The non- 
weighted fit for afl the data at z > lis (0.76+0. 18)x(l-Hz)2°"^*'-25. 

Interestingly, some earlier numerical simulations and theo- 
ries with a quasar-only UV background have shown that there 
should be a break in the dn/dz evolution at z ~ 2 due to the de- 
crease in the quasa r number density, thus less available H i ionis- 
ing photons ( Theuns et al.lll998l;lDave et al.lll999HBianchi et al.l 
2001). The green dot-dashed curve in Fig. |4 shows one of 
such predicted dn/dz evolutions bv iDave et al. (Il999l) . which 
outlines the Weymann et al. dn/dz reas o nably well. However, 
more recent simulations by iDave et aP (|2010) predict differ- 
ent dn/dz evolutions. These simulations are based on the vari- 
ous galactic wind models and the UV background contributed 
both by quasars and galaxies. The red dotted and the blue 
dot-dot-dot-dashed curves at z < 2 illustrate their predicted 
dn/dz based on momentum-driven wind and no- wind models, 
respectively. These newer simulations predict that dn/dz con- 
tinuously decreases with decreasing redshift. Their momentum- 
driven wind model agrees reasonably well with the observations 
by HSTfSTlS with the H i absorber selection imposed (measure- 
ment errors less than 40% and b < 40 km s '), but not with the 
Weymann et al. data. A better, uniform dataset from HSTfCOS 
observations should resolve this discrepancy at z < 0.5. 

For the column density interval for stronger absorbers 
logA^Hi = [14.5, 17.0], our data shows that the evolution con- 
tinues to follow the empirical power-law with y = 3.02 + 0.48 
(see Table |2]l. However, the scatter between different sightlines 
is large as stronger absorbers are rare at all redshifts (Dave e t aP 
2010). There are more than 3cr difference between the lowest 
dn/dz sightline and the highest dn/dz sightline atz~2. lKim et alJ 
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Fig. 5. The number density evolution of the Lya forest in the 
column density range logA^Hi = [13.1, 14.0]. All the symbols 
have the same meaning as in Fig.|4] 



( l2002h discuss the possibility on whether the column density 
evolution flattens out at z < 2.5 (log(l + z) < 0.55) for this col- 
umn density interval. Even though more data points are available 
in this study, this question cannot be conclusively answered and 
more data covering lower redshifts are required. 

The line number density evolution for low column density 
systems in the range of logA^Hi = [13.1, 14.0] is presented in 
Fig- IS] Similar to Fig.|4] it suggests that the flattening of dn/dz 
at z < 1.5 might continue at the lower column density range. 
However, the sightline variation at z < 0.4 is larger at this col- 
umn density range. This is in part caused by different analysis 
methods and diff'erent S/N STIS data used by different studies. 
For example, the number density measured in the STIS spec- 
trum toward PKS0405-123 is different between the Williger et 
al. (2006) work (filled purple upside-down triangles) and the 
Lehner et al. (2007) work (two of open red upside-down trian- 
gles). Again the results from our data agree well with previous 
results found in the literature at z > 1.5. The linear regression 
to our data at z > 1 .5 gives y ^ 1 .67, comparable to the fit in- 
cluding all available literature data points at z > 1 .0. However, 
these resu lts do not c ompare well with the linear regression ob- 
tained bv iKim et al] (|2002) with y = 1.18 + 0.14 (the dotted 
line), a shallower dn/dz evolution. This discrepancy arises due 
to their rather small sample size at z < 2.5 and more severe line 
blending at higher redshifts. Given a larger cosmic variance at 
lower redshifts, the sample size becomes more important. At the 
same time, line blending at high redshifts makes the detection 
of weak absorbers difficult. This incompleteness effect has been 
shown to underestimate the line number density of low column 
density systems at logA^Hi - [13.3, 13. 6] by ~ 17% at z ~ 3 
(log(l +z)~ 0.6) and by ~ 35% at z ~ 4 (iGiallongo etanil996l). 
Both effects tend to flatten the evolution observationally from 
its true value. In addition, the robust estimate of the exponent y 
requires a large z leverage. 

Even though there are not many sightlines covering 0.5 < 
Z < 1.5, we calculated the mean dn/dz from all the combined 
Hi fitted line lists including the literature data in Fig. |6] This 
mean dn/dz is not an averaged value of the individual sightlines. 
The literature data used in the combined line list include all the 
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Fig. 6. The mean number density evolution of the Lyo' forest. 
The vertical error bars give the Icr Poisson error, while the x- 
axis error bars show the redshift range covered by each data 
point. The data point with the dotted error bar indicate the 
iJanknecht et al.l (l2006h data. The dashed straight lines mark the 
fit to the mean data excluding the .Janknecht et aL (.2006.) data. 



quasar sightlines shown in Figs. |4]and|5] except the HST/FOS 
Weymann et al. (1998) data, the WilHger et al. (2006) data and 
the Savaglio et al. (1999) data. The HST/FOS data was excluded 
since they were based on the equivalent width measurements, 
while the Williger et al. (2006) data suffered from noise features. 
The Savaglio et al. (1999) result is from a single sightline and 
provides the only data point besides the Janknecht et al. (2006) 
data at z ~ 1. Although the Janknecht et al. (2006) data also 
suffer from noise, they were from 9 sightlines. We opted to use a 
result based on the analysis of multiple sightlines from a single 
study. This helps to reduce any systematics caused by combining 
results from different studies at z ~ 1 . For z < 0.4, the systematic 
uncertainty is larger since the line lists used are produced by 
different studies. 

At logA^Hi = [13.1, 14.0], there might occur a flattening at 
z~ 1, if the Janknecht data were included. At log A^hi - [14, 17], 
a single power law with y - 1.61 ± 0.12 does not give a good 
fit at < z < 4, regardless of the inclusion of the Janknecht et al. 
data. It remains to be seen whether a single power law fits the 
dn/dz evolution for both high and low column density ranges 
at z = 0.4. It should be noted that the dn/dz of Lyman limit 
systems with a column density of logA^Hi - [17.2, 19.0] does 
not fit to a single power law. It shows a slower evolut ion at z < 2 
and evolves rapidly at z > 2 ('Procha ska et al.ll2010l) . while the 
dn/dz ofdampedLya systems with log A^hi = [20.3,22.0] shows 
a single pow er-law evolutio n with a slope y = 1 .27 ± 0. 1 1 at 
< z < 4.5 ( Rao et al.ll2006h . 

Our results indicate that higher column density forest sys- 
tems evolve more rapidly than low column density systems and 
the number density of high column density systems decreases 
faster with decreasing redshift. The increase in the scatter at red- 
shifts z<2.5 might indicate the transition point where the evolv- 
ing number density changes into a non-ev olving one, as is pre- 
dict ed in earlier nume rical simulations by iTheuns et al.l (Il998t) 
and lDave et alJ(ll999l) . 
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3.2. The differential column density distribution function 

The differential column density distribution function (CDDF) is 
defined as the number of absorbers per unit absorption distance 
X{z) and per unit column density A^hi- The absorption distance is 
calculated using Eq. [1] Empirically, the differential distribution 
function is reasonably well described by a single power law at 
z ~ 3 atlogA^Hi = [13,22] as 



1 + 6 



An 



An 



Ni 



Hi' 



(4) 



where (dn/(dA^Hi dX))^ gives the normalisation point of the dis- 
tribution function and /3 denotes its slope. However, the detailed 
shape of the differential column density distrib ution function is 
dependent on the A^hi column density range dProchaska et al.l 
120101: lAltav et allbOl 1,) . It shows a flattening around the tran- 
sition from the forest to the Lyman limit systems at A^hi at 
log A^Hi ~ 17. Then it shows a steepening at log A^hi ~ 20 where 
a transition occurs from the sub-damped Lyo- systems to the 
damped Lya systems. 

In Fig. |7] we present the results using the Lya-only fits at 
1.9 < z < 3.2. Note that the redshift range used for the CDDF 
analysis is different from the one used for the dn/dz analysis in 
Section [TTI The total absorption distance X(z) at 1.9 < z < 3.2 
is 21.8165. The binsize of logA^Hi - 0.25 is used at logA^Hi - 
[12.0, 15.0], then tiie binsize of 0.5 at logA^Hi = [15.0, 18.0]. 
To in crease the column densi ty co verage, we include re sults 
from iNoterdaeme et al.1 (l2009l) and 'O'M eara et all (|2007l) for 
logA^Hi > 19 at z ~ 3 from flie SDSS II DR7 datfl The top 
X-axis is in units of th e gas overdens ity 6 which was computed 
according to Eq. 10 bv lSchav3(l2001b 



2.7 X 10'^ (1 +6)^'^ X 



7.-0.26 T^-l 
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Here, the gas temperature T is assumed to be T = T4 x 10'^ K, 
the photoionisation rate F - Fn x 10"^^ sec"'. The parameter 
fg denotes the fraction of mass in gas. The IGM gas tempera- 
ture is assumed to be governed by the effective equation of state 
T - Tp jl + df'^, where T p is the temperature at the cosmic 
density dHui & Gnedinlll997h. For Fi? and 7, we interpolated re- 
sults obtained by Bolton et al] (l2008h . We assumed th at Tq is 
2 X lO'* K, fg = 0.16, and Q^/z^ = 0.0227 (ISchavell20 0l'). As the 
same overdensity corresponds to a different A'hi at different z, 
the overdensity plotted in Fig.Qis at the mean redshift, z = 2.55. 

W e compa re ou r results wit h the observations by 
IPetitiean et all ( Il993l) and Hu et alj (Il995l) . Our results are 
in good agreement with the IPetitiean et al.l (1 19931) data over 
the whole column density range down to logA^Hi ~ 13.5, 
following a power law at logA^Hi = [12.75,14]. At smaller 
column densities logA^Hi < 12.75, the CDDF starts to devi- 
ate from a power law due to the sample incompleteness for 
weak absorbers (iKim et alJll99'7i) . From the linear regression, 
we find log(dn/(dA^HidX))o = 7.34 ± 0.42 and a slope of 
/3 = -1.43 + 0.03 for the logA^Hi = [12.75, 14] range (the solid 



^ The plotted d ata points from the l iteratur e are the reported ones in 
each study. Both O'Mea ra etal] (12007) and Noterdaeme et al.| ( l2009h 
used the same cosmology as ours, while Petitiean et al. ( 1993) used the 
qo = cosmology. The absorption distance X(z) in our cosmology is 
about 6% smaller than theirs. Since the CDDF uses the logarithm value 
of X(z), the difference in the CDDF is negligible even without convert- 
ing their CDDF to our cosmology. 
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Fig. 7. The differential column density distribution at 1 .9 < z < 
3.2 using the Lyc-only fits. Both black (logA^Hi ^ 12.75) and 
grey (logA^Hi < 12.75) data points show the results from our 
quasar sample. The grey data points mark the column densi- 
ties that are affecte d by incompleteness. The stars are the data 
points obtained by IPetitiean et alj (Il993h . The filled circles at 
logA^Hi > 20.1 and at 19 < logA^Hi < 20.1 are from the SDSS 
II DR7 at < z >= 3 .02 bv INoterdaeme et al.l (|2009|) and from 
lO'Meara etal] (l2007h at < z > = 3.1, respectively. The solid fine 
gives the power law fit to our data for lo g A^ht = [12.75, 1 4]. The 
dotted line represent the fit obtained bv lHu et al] (1 19951) . while 
the dashed hne re presents a theoretical prediction at z ~ 3 by 
lAltavetalJ (1201 Ih . The vertical error bars indicate 1 cr Poisson 
errors, while the x-axis error bars show the A^h i range covered 
by each data point. Gas overdens i ties on the top x-axis are com- 
puted using Eq. 10 from lSchavd (12001 1) at z = 2.55 (see text for 
details). 



line). Thi s result i s sligh tly lower than /3 = -1.46 (no errors 
given) b vlHu et al.l(ll995h (flie dotted line) oryS = -1 .49 ± 0.02 



bv iPetitiean et al 



The distribution function becomes steeper at logA^Hi > 15, 
then be comes shallower at higher logjVHi, a s previously ob- 
served (IPetitiean et all 119931: iKim et alJ 119971: IProchaska eTaH 
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This result agrees well with the theore tical prediction at 
Z ~ 3 (the dashed line) bv lAltav et d!k(l201 ih at logA^Hi Sol6, 
but starts to show a noticeable disagreement at the 1-3 cr level 
at logA^Hi = [16, 18], in part due to the lack of enough high- 
column density systems in our small sample. We will address 
the shape of the CDDF in more detail in the next section using 
results from the high-order fit sample. 

4. Analysis using liigher-order Lyman lines 

In the last section, we checked the Lya absorber number den- 
sity evolution and the differential column density distribution 
obtained from the Lya-only fits for consistency with previous 
studies. The analysis is now revisited with the results from the 
Voigt profile analysis including the higher order transitions at 
1 .9 < z < 3 .2, hence a sample with a more reliable A^h i ■ Therefore, 
it can be established whether the dip seen in the differential col- 
umn density distribution at logA^Hi between 14.5 and 18 is a 
physical feature or just an imprint of uncertainties in A^hi- All 
the results from this section are tabulated in Appendix A. 

4.1. The mean number density evolution 

We now revisit the line number density evolution using the high- 
order Lyman sample, as described in the previous section. On a 
quasar by quasar analysis we determine dn/dz for a low column 
density range of logA^Hi - [12.75, 14.0] and for high column 
densities of logA^Hi = [14, 17]. The lower column density range 
is chosen in such a way that the part of the differential column 
density distribution function which follows a power-law is cov- 
ered, whereas the logA^Hi - [14, 17] interval covers those sys- 
tems responsible for the dip in the column density distribution 
function. 

The results are presented in Fig. |8] Linear regressions from 
the data are obtained and the resulting parameters are sum- 
marised in Table |2] Similar to the previous analysis, the line 
number density shows a decrease with decreasing redshift. No 
significant differences between the two different fits are present, 
even though the total redshift coverage used for the high-order 
fit is about 20% smaller In the case of the logA^Hi - [14, 17] 
interval, the slope of the power law steepens from the Lyor-only 
slope of 7 = 3.40 + 0.36 to y = 4.91 + 0.53 for the high-order 
fit. This is in part caused by that the number of high column den- 
sity absorbers is larger in the high-order fit sample. However, the 
slopes of the two samples are still in the 2cr uncertainty range, 
rendering the two results consistent to each other Similar results 
are obtained for the log A^hi - [12.75, 14.0] range. The slope for 
the high-order fit increases from 7 = 1.13+0.16fortheLyQ'-only 
fit to 7 = 1.38 + 0.22. Again, the results from the two samples 
agree within the Icr uncertainty range. 

In previous studies the number density evolution has been 
usually derived on a quasar by quasar analysis. Previous studies 
did not have enough quasar sight lines available to sample the 
number density evolution at smaller redshift interval Az, without 
suffering from small number statistics. Our sample of 18 high- 
redshift quasars is characterised by a large redshift distance cov- 
erage in the redshift range of 1.9 < z < 3.2 (see Fig.|2]i. As a 
result, a large number of absorption lines is available for small 
redshift intervals to combine the individual quasar line lists into 
one big sample. However, due to the larger cosmic variance at 
low redshifts from the structure formation, the redshift bin size 
should not be too small. From this combined sample, the evo- 
lution of the mean number density is derived in redshift bins of 
Az - 0.26, starting fromz - 1.90. 




4.0 
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Fig. 8. The line number density evolution derived on a quasar 
by quasar analysis using the high-order Lyman sample for col- 
umn density intervals of log A^hi - [12.75, 14.0] and logA^Hi = 
[14, 17]. The vertical error bars indicate 1 cr Poisson errors, while 
the X-axis error bars show the redshift range covered by each 
sightline. The straight lines denote results from a linear regres- 
sion to the data with parameters given in Table |2] The data are 
tabulated in Table lA. II 



Results of the combined line number density evolution are 
shown in Fig. |9] for identical column density ranges as used in 
the quasar by quasar analysis. Error bars have been determined 
using the bootstrap technique. For comparison, results using the 
LycK-only fits are overplotted as grey open circles for the high 
column density bin. 

The high column density results are similar to the ones 
obtained from the Lya-only fits. The number density itself is 
higher in the high-order fits, since some strongly saturated sys- 
tems break up into multiple, strong components in the high- 
order Lyman transition. In addition, three absorbers (two toward 
HE0940-1050 and one toward Q0420-388) were found to be 
a Lyman limit system with logA^ei > 17 in the Lyc-only fit. 
Therefore, these systems were not included in the Lya-only re- 
sults. However, these Lyman limit systems break up into mul- 
tiple weaker components in the high-order fit and contribute to 
the number count in the high order fit analysis. However, the dif- 
ferences between the two samples are smaller than the statistical 
uncertainties. 

At low column densities, no noticeable differences between 
the two samples are observed, as expected. 

Again, linear regressions have been determined and their pa- 
rameters are given in Table |2] At logA^Hi = [12.75, 14.0], the 
slope of our combined sample is 7 = 1.28 + 0.24, similar to 

7= 1.38±0.22fromthequasarby quasar analysis. At log A^H I = 
[14, 17], the slope of the combined sample 7 - 4.65+0.66 is also 
similar to 7 = 4.91 + 0.53 obtained from the quasar by quasar 
analysis. The slopes from both analyses of our high-order fit 
sample at logA^Hi = [14, 17] are steeper than the ones obtained 
from the Lya-only fit sample. In particular, the ones from the 
combined sample differ more than 3cr. This difference is mainly 
caused by that the redshift range used for the combined sam- 
ple is different for two analyses. For the Lycc-only fit, the mean 
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Fig. 9. The mean line number density evolution of the com- 
bined sample as a function of redshift using the high-order 
Lyman-series sample for column density intervals of log A^hi = 
[12.75, 14.0] and logA^Hi = [14, 17]. The sample is binned in 
redshift with Az - 0.26, starting from z = 1.90. The vertical er- 
ror bars indicate Icr Poisson errors, while the x-axis error bars 
show the redshift range covered by each data point. For compari- 
son, the results of the Lyc-only fits (grey open circles) are shown 
for the log A^Hi = [14, 17] interval. The straight solid lines denote 
results from a linear regression to the binned data. Two dashed 
lines represent the mean number density evolution of the Lya- 
only fit sample for logA^ei - [14, 17] (logdn/dz - (-0.41 + 
0.35)-H(4.14 + 0.63)xlog(l -Hz))andforlogA?Hi = [12.75,14] 
(logdn/dz = (1.96 + 0.13) + (1.12 + 0.24) x log(l + z)), respec- 
tively. Exactly same redshift range was used for both fit samples. 
The data are tabulated in Table IA.2I The parameters of the fits 
are given in Table |2] 



d«/dz is derived for < z < 4, while for the high-order fit it is 
restricted to 1.9 < z < 3.2. 



4.2. The differential column density distribution function 

Using the high-order fits, we have derived the differential col- 
umn density distribution function (CDDF) for 1 .9 <z < 3.2, anal- 
ogous to Section [J!2] In Fig.[TO]we show the results for the entire 
redshift range. As in Fig.Q the binsize of log A^hi =0.25 is used 
at logA^Hi = [12.0, 15.0], then the binsize of 0.5 at logA^Hi - 
[15.0, 18.0]. The total absorption distance is X{z) = 21.8165, the 
same value used for the Lya-only fit CDDF analysis. As with 
the Ly g-on ly fits, we inc luded observations bv lNoterdaeme et aO 
(I2009h and lO'Meara eTal., (,2007) . 

The high-order fit results show a power law relation which is 
almost identical to the results of the Lyc-only fits. As with the 
Lya-only fits, the differential column density distribution func- 
tion shows a deviation from the empirical power law at column 
densities between 14 < logA^Hi < 19. Since the column density 
distribution deviates from a single power law at logA^Hi ~ 14, 
we have individually fitted power laws to four column den- 
sity intervals of logA^Hi = [12.75, 14.0], [14, 15], [15, 18], and 
[12.75, 18.0] at 1.9 < z < 3.2, characterising the shape of the dis- 
tribution function. The resulting parameters are listed in Table[3] 
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Fig. 10. The differential column density distribution at 1 .9 < 
z < 3.2 of our quasar sample using the high-order Lyman fit. 
Black and grey data points show the results from our quasar 
sample. The grey data points below logA^Hi < 12.75 mark the 
column densities that are affected by incompleteness. The grey 
data points above log A^hi > 12.75 represent the results from the 
Lyof-only fit. The dashed lin e represents a theoretical prediction 
at z ~ 3 bv lAltav et al.l (1201 lb . The vertical error bars indicate 1 cr 
Poisson errors, while the x-axis error bars show the A^hi range 
covered by each data point. The filled circles at log A^hi > 20.1 
and at 19.0 < logA^Hi < 20.1 a r e from the SDS S II DR7 at < 
z>= 3 .02 bv lNoterdaeme etal.1 (|2009() and from lO'Meara et al.l 
(I2OO7I) at <z> = 3.1, respectively. The solid line gives the power 
law fit to our data for logA^Hi - [12.75, 14.0]. The overdensity 
plotted on the top x-axis is calculated at z = 2.55. 



AtlogA^Hi = [12.75, 14.0], the linearregression yields a nor- 
malisation point of log(dA^/(dA^HidX))o - 7.41 + 0.42 and a 
slope of p - -1.44 + 0.03. This result is almost identical to 
the Lya-only fit, since differences between the LyfT-only and the 
high-order fits start to be significant at logA^Hi > 15 (see Fig. 
[Tol l. The high-order fits show a larger number of absorbers at 
15 < logA^Hi < 17 than the Lya-only fits. However, at higher 
column densities, the number of absorbers is lower for the high- 
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order fits than for the Lya-only fits. This again indicates the 
breaking up of high column density systems into multiple lower- 
A^Hi ones when including higher transitions than Lya. For the 
entire redshift sample, the slope becomes steeper from — 1.44 
to ~ -1 .67 at log A^Hi = [14, 15]. Then at the higher column den- 
sity range logA^Hi = [15, 18], the slope becomes shallower to 
— 1.55, a trend shown in the numerical simulation (the dashed 
line) by Altay et al. (2011) in Fig.fTOl 

In order to determine the redshift evolution of the differential 
column density distribution, we split the sample into two redshift 
bins: z = [1.9,2.4] and [2.4,3.2]. Fig. [lU shows the CDDF at 
two different redshift bins, where we overplot the power-law fit 
at log A^Hi - [12.75, 14.0] for each redshift bin as the solid line. 
We also overplot the results of the power-law fit to the entire 
redshift range 1.9 < z < 3.2 at the same column density range 
as the dashed line. For the redshift intervals of z = [1.9,2.4] 
and [2.4, 3.2], the absorption distance is X(z) - 11.8049 and 
10.01 16, respectively. 

Unfortunately, the uncertainties in the power-law fit param- 
eters at each redshift bin do not allow us to constrain the shape 
of the distribution as a function of redshift reliably. Comparing 
the slope of the linear relations shows that the CDDF becomes 
slightly steeper at low redshift for log A^hi - [12.75, 14.0], from 
y = -1.38 + 0.04 at high z to y = -1.49 + 0.04 at low z. 
However, the slopes are still consistent within 2cr, i.e. no signif- 
icant CDDF evolution, cf. ^Villiger et al. (2010). They are also 
consistent with the result from the entire redshift range within 
Icr. 

Let us now focus on column densities above logA^ei > 14.0. 
From Fig. [To] we have seen that the diff'erential column density 
distribution deviates from the power law form for column den- 
sities logA^Hi > 14.0. The lower panels of Fig. [TT] show the dif- 
ference between the observed CDDF and the power-law fit to 
the CDDF for the entire redshift range (the dashed lines). The 
entire redshift fit was used since the comparison requires an ab- 
solute reference. From the lower panels, it is clear that the de- 
viation from the power law is stronger for the low redshift bin. 
At the same time, the deviation column density above which the 
deviation starts to be noticeable is lower at low redshift, from 
log A^Hi ~ 14.7 at 2.4<z<3.2 to logA^Hi- 14.0 at 1.9<z<2.4. 

Note that no such break in the CDDF has been seen in the 
log A^Hi = [12.5, 16] at z < 2, cf. Fig. 5 of IWiUmer et al.l (1201 Oh 
at z ~ 0.08 and Fig. 9 of 'Ribaudo et alj (l201l!) at z < 2. Both 
works also found a steeper CDDF slope of ~1.75. Some of the 
discrepancy is caused by the different fitting methods, the H i se- 
lection criterion discussed in Section lTTI and the column density 
range over which the powe r law was performed. On the other 
hand,'Prochaska et al. (2010) found a more significant dip in the 
column density distribution function at logA^ei = [14,19] at 
z ~ 3.7 (similar to the Altay simulation at z ~ 3 indicated by 
the dashed curve in Figs. iTl and [TOll. However, the dip shown at 
z ~ 3 in Fig. [TT](the high redshift bin) is not as strong as the one 
predicted by the Altay simulation, although both results are still 
considered to be consistent within 2cr. These diff'erences could 
be simply due to our small sample size, or due to the different 
analysis method or due to the strong CDDF evolution between 
z~4 and z~0. 

Note that the dip shown in Fig. [TT] is not caused by self- 
shielding. Self-shielding causes the number density of absorbers 
to increase. Self-shielding becomes important at logA^Hi ^ 16 
and its eff'ect becomes evident at logA^Hi ^ 17 with a shal- 
lower slope than the extrapolated one at the lower logA^ei 
dAltav et al.l 1201 ll) . However, the dip in discussion occurs at 
logA^Hi = [14.5, 17.0] compared to the extrapolated power-law 



slope at logA^Hi - [12.75, 14.5]. In addition, the deviation A^hi 
from this single power law starts at log A^hi ~ 14.5, where self- 
shielding has no effect. 

5. Characteristics of the metal enriched forest 

The discovery of metal lines which are asso ciated with H i ab- 
sorber in the Lya f orest, such as C iv or O vi dCowie et al.ll 19951: 
ISongailai 1998; Sc have et al.l l2000ah . have raised the question 
of how the IGM has been metal enriched. As the forest has 
a high temperature and a low gas density, it is not likely to 
form stars in-situ. Metals should be transferred from galax- 
ies by e.g. galactic outflows (Agum'e et al. 20()Tt ISchave et al.l 
l2003t lOppenheimer & Davell20()6l)T In recent years, studies on 
galaxy-galaxy pairs at high redshift have revealed some evidence 
that metals associa ted with the Lya forest r eside in the c i rcum- 
galactic medium (lAdelberger et al.l 120051: ISteidel et all 120101: 
Rudie et aL,2012) . In this interpretation, the metal-enriched for- 
est cannot be called the IGM in the conventional sense and is 
likely to show a different evolutionary behaviour compared to 
the metal-free forest. In order to learn more about these enriched 
hydrogen absorbers, we characterise C iv enriched H i absorbers 
in this section by determining their number density evolution and 
differential column density distribution. Note that we excluded 
Q0055-269 and J2233-606 for both the C iv enriched forest and 
the unenriched forest samples in this section, as their C iv region 
has a much lower S/N of ~ 40 per pixel compared to the other 16 
quasar spectra whose S/N is greater than 100 per pixel in most 
C IV regions. Due to the wavelength gap caused by the UVES 
dichroic setup, the C iv redshift coverage is shorter than the H i 
coverage for Q0420-388, HE0940-1050 and HE2347-4342. 
We excluded the +200 km s"' region from the wavelength gap 
and included the C iv region only when it covered both doublets. 
The excluded C iv redshift range for these three quasars is listed 
in Table [T] In this section, we used the column density and b pa- 
rameter of H I from the high-order Lyman fit, unless stated other- 
wise. All the results from this section are tabulated in Appendix 
A. 



5.1. Method 

Unfortunately there is no one-to-one relation between H i lines 
and C IV lines. Fig. [12] shows a velocity plot (the relative ve- 
locity centered at the redshift of an absorber vs normalised 
flux) of a typical C iv-enriched H i absorber in the spectrum of 
HEl 122-1648. The vertical dashed lines indicate the velocity 
of individual Hi components. Not all Hi lines can be directly 
assigned to one or only one C iv component. For example, the 
H I component at -36.02 km s"' could be associated either with 
the first Civ component at -51.72 km s ' or with the second 
one at -19.97 km s or with both. A general trend is that the 
associated C iv features show an increased number of velocity 
components as A^hi increases. The absorption line centers of H i 
and C IV lines often show velocity differences as well, indicating 
that the H i-absorbing gas might not be co-spatial with the C iv- 
producing gas. Therefore, we apply a simple assigning method 
to our fitted absorber line lists, in order to determine if an H i 
absorption line is associated with C iv. 

We consider an H i absorber to be metal enriched if a C iv 
line with A^c iv greater than a threshold value exists within the 
velocity range +Avciv centered at each identified Hi line. The 
threshold A^civ should be large enough not to be affected by the 
incompleteness of weak C iv detection, but not too large so that 
there are enough C iv enriched absorbers to have a meaningful 
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Table 3. Linear regression results for the differential column density distribution as a function of redshift and column density using 
the high-order fit. Here the normalisation point log (dn/(dA^Hi <iX)\) is denoted by B. 
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Fig. 11. Upper panels: The differential column density distribution as a function of redshift. Both black and grey data points show 
the results from the high-order Lyman sample. The grey data points at log A^hi < 12.75 mark the column densities that are affected 
by incompleteness. The vertical error bars indicate Icr Poisson errors, while the x-axis error bars show the A^hi range covered by 
each data point. The black solid line gives the power law fit for log A^hi = [12.75, 14.0] at each redshift bin, whereas the dashed line 
is the fit to the z = [1.9, 3.2] redshift range (see Fig. [TOt . The overdensity plotted on the top x-axis is calculated at the mean z for 
each redshift bin. Lower panels: the residuals from the power law fit from the entire redshift range at log A^hi = [12.75, 14.0]. 



statistics. This method can assign one Hi component with mul- 
tiple C IV components and vice versa. As we are not concerned 
with the one-to-one relation between A^c iv and A^h i of each H i 
component, but the existence of the C iv line for a given search 
velocity range, the multiple assigning of the same component 
does not affect the results. 

Two arbitrary choices of Avciv are considered: a conserva- 
tive narrow range of +10 km s"' (a minimum b value of a single 
Lya absorption line is roughly 20 km s ') and a more generous 
interval of + 100 km s"' . 



Fig. [T3] shows the Civ column density distribution function 
at 1.9 < z < 3.2 from our sample (black filled circles). For 
comparison, other results from the literature are als o included: 
red fil led triangles and blue filled diamonds from Pichon et al.l 
(I2003h at 1.5 < z < 2.3 and 2 . 3 < z < 3.1, respectively, and green 
open squares from lSongailal (1200 ll) at 2.90 < z < 3.54. The turn- 
over seen in green open squares is due to the incompleteness 
effect, i.e. not all weak C iv can be detected due to noise. 

Similar to the Hi density distribution, the Civ CDDF 
does not fit with a single power law over a large A^civ 
range. The Pichon et al. result even suggests that the Civ 
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Fig. 12. Example of a velocity plot (a relative velocity vs 
normalised absorption profile plot) of Hi and associated Civ 
detected in the z - 2.352743 absorber in the spectrum of 
HEl 122-1648. The strongest C iv component is set to be at the 
zero velocity. The observed spectra are plotted as a histogram, 
while Voigt-profile fits are as a smooth curve. Thick red curves 
are the combined fit from individual components. The heavy tick 
marks above the profiles indicate the velocity centroid of each 
component. Non-negligible blends by other ions are indicated in 
gray. The h value (in km s ') and log A^hi with the VPFIT fitting 
errors are displayed next a tick mark indicating the center of the 
component. 



density distribution might have a non-linear functional form. 
At logA^civ = [12.25,13.5], a single power-law fit gives 
logdA^/(dA?civdX) = (6.60 + 1.23) h- (-1.47 + 0.10) x logA^civ 
(the solid line). At logCiV - [12.25, 15.5], a single power law 

is log dA^/(dA^c ivdX) = ( 1 1 .4 1 + 1 .6 1 ) -H (- 1 .85 + 0. 1 3) X log A^c IV 
(the dotted line). If the solid line is taken as a reasonable CDDF 
since it fits the low-A^c iv CDDF better, our C iv detection can be 
considered complete at log A^civ ^ 12.2. 

Another way to look at whether our A^c iv completeness limit 
is reasonable is with the column density-/? value diagram. As 
seen in the 7th column of Table [1] the S/N differs for differ- 
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Fig. 13. The C iv column density distribution. Filled circles are 
our results at the redshift range used for the high-order fit C iv- 
enriched Hi sample at 1.9 < z < 3.2. The CCD gap in the Civ 
region was account ed for Red filled tr iangles and blue filled 
diamonds are from iPichon et al.l (l2003l) at 1.5 < z < 2.3 and 
2.3 < z < 3.1, r espectively. Green open squares are taken from 
ISongailal(l2001h at 2.90 < z < 3.54. The vertical error bars indi- 
cate Icr Poisson errors, while the x-axis error bars show the A^hi 
range covered by each data point. The black dotted line shows 
the linear regression to filled circles at log Ncw - [12.25, 15.5] : 
logd«civ/(dA?civdX) = (11.41 + 1.61)+(-1.85+0.13)xlogA^civ 
The solid line is the power law fit at logA^civ - [12.25, 13.5] : 
logd«civ/(dA?civdX) = (6. 60 +1.23) -H(- 1.47 + 0.10) X log A^c IV 
The turn-over at logA^civ ~ 12.5 shown in green data is simply 
due to the incompleteness for weak C iv. Similarly, the turn-over 
seen at log A^civ ~ 12.1 in our data is also due to the incomplete- 
ness. 



ent sightlines, and changes even along a single spectrum. This 
makes it extremely difficult to quantify the correct 3cr detection 
limit for a dataset containing spectra with different S/N. 

Fig.ll4lshows the logA^civ-^civ diagram at the two redshift 
bins. The vertical heavy dot-dashed lines mark log A^civ = 12.2. 
In the upper panel, two heavy dashed lines show a 3o- detection 
limit for a spectrum with S/N = 120 (the left side) and 90 (the 
right side, an approximate lowest S/N) per pixel, respectively. In 
the lower panel, the heavy dashed line is a 3cr detection limit 
for S/N - 120. Absorption lines at the left-side of the detection 
limit, i.e. broader and weaker lines, can be only detected for S/N 
greater than the given S/N. Overlaid as a histogram is the dis- 
tribution of the number of Civ lines with logA^civ ^ 12.2 as a 
function of hcw- For the distribution, the zero base is set to be 
log A^c IV = 12.2. Thick ticks above the distribution mark the me- 
dian bcw- There is no correlation between A^civ and /^civ above 
the S/N detection limit at all of the reasonable expected hcw 
values. 

At 2.4 < z < 3.2 (the upper panel), the 3cr bcn detection 
limit is 23.6 (13.4) km s"' for S/N = 120 (90) at log A^civ ~ 12.2. 
The total wavelength coverage of C iv at the high redshift bin 
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Fig. 14. Line width vs. column density for the Civ absorption 
lines along 16 sightlines excluding Q0055-269 and J2233-606 
at the two redshift bins. Error bars are fitting errors from the 
VPFIT profile fitting. In the upper panel, two heavy dashed lines 
delineate a 3cr detection limit for a spectrum with S/N = 120 
and S/N - 90 per pixel. At 2.4 < z < 3.2, most spectra show S/N 
greater than 90. In the lower panel, the heavy dashed line shows 
a 3cr detection limit for S/N - 120. Broader and weaker absorp- 
tion lines at the left of the detection limit are missed in lower 
S/N spectra. The vertical dotted line indicates the adopted low 
A^civ bound of log A^civ - 12.2 above which the incompleteness 
does not affect the Civ detection significantly. The histogram 
shown with the base at logA^civ = 12.2 is the number of Civ 
lines as a function of bcw with the bcw binsize of 2 km s 
Thick ticks above the number distribution mark the median bcn- 
for logA^civ ^ 12.2. The total number of Civ lines is 194 and 
171 at 2.4 < z < 3.2 and 1.9 < z < 2.4, respectively. Among 
them, 138 and 122 lines have logA^civ ^ 12.2 at the same red- 
shift range. 



is ~ 3192 A. For about half of the spectra there is contami- 
nation from weak telluric lines in < 10% of the Civ region. 
This contamination prevents isolated weak C iv lines from be- 
ing detected, however, can be treated as a lower-S/N region. 
Including the telluric-contaminated region, the wavelength cov- 
erage with S/N < 120 is about 1018 A. In the Civ wavelength 
region with S/N > 120, the total number of Civ lines with 
logA^Civ - [12.2, 12.3] is 8. Out of those 8, none has bciv>23.6 
km s^' . It is possible that a large fraction of C iv has a ^civ value 
greater than 23.6 km s and therefore, would be completely 
missed even in the high-S/N spectra analysed here. However, as 
clearly seen in the upper panel of Fig. [14] the bciv distribution 
at logA^civ > 12.2 shows that only 9% of Civ has bcn > 23.6 



km s ' . If a large fraction of C iv lines were broader regardless 
of A^civ, the region around log A^civ ~ 12.4 and bcn ~25 km s"' 
in Fig. [T4| should have been more crowded. Therefore, it is not 
likely that many weak Civ lines with bciv > 23.6 km s"' have 
been missed for S/N > 120. 

Only 2 out of 8 have ^civ ^ 13.4 km s"' at logA^civ = 
[12.2, 12.3]. In other words, these 2 C iv lines would have been 
missed in the S/N < 120 region. One is a single isolated line, 
while the other is part of a multi-component C iv complex. We 
assumed that the number of C iv lines with log A^civ ~ 12.2 and 
bciv ^ 13.4 km s"' is 2 in the wavelength range of 2174 A, 
i.e. the total wavelength range with S/N > 120. If we assume 
that weak Civ lines have a negligible clustering, about 1 (or 
2 X 1018/(3192 - 1018) = 0.9) Civ line with logA^civ ~ 12.2 
and bciv > 13.4 km s"' could have been missed in the C iv forest 
region with S/N <120. 

A total of 5 Hi lines with logA^Hi = [12.8, 17.5] is found 
within +100 km s ' centered at these two Civ lines. The total 
number of high-order-fit H i lines in the H i forest region corre- 
sponding to the S/N > 120 Civ forest region is [265, 363, 233, 
120, 50, 27, 5] for logA^Hi = [12.75-13.00, 13.0-13.5, 13.5- 
14.0, 14.0-14.5, 14.5-15.0, 15.0-16.0, 16.0-17.0], respectively 
Among them, a negligible number of H i lines, [0, 0, 2, 0, 1, 0, 
1], is associated with these two C iv lines for the same A^hi range, 
or less than 2%. The remaining one Hi line has logA^Hi ^ 17.0 
as the associated C iv line belongs to a C iv complex of a partial 
Lyman limit system. Although the number of undetected weak 
and broad C iv lines in the S/N < 120 region is a very rough esti- 
mate, less than 2% of the Hi lines would be mis-classified as the 
unenriched forest due to the incompleteness at log A^civ ~ 12.2. 

The situation becomes more complicated in the low redshift 
bin, where the variance of the S/N limits of individual spectra 
is much higher than in the high redshift bin. If a similar logic 
were applied to, the total C iv coverage is 5485 A, and the one 
with S/N < 120 is 2719 A. In the S/N > 120 Civ region, there 
is a total of 10 Civ lines with logA^civ = [12.2, 12.3]. Out of 
10, 6 lines have bciv S 13.4 km s"', the maximum bcw value 
to be detected for a line with logA^civ - [12.2, 12.3] in a S/N 
- 90 spectrum. Among those 6 C iv lines, two C iv lines are part 
of a two-isolated-component complex, with the rest being part 
of a multi-component complex. Since stronger H i lines tend to 
be associated with a C iv complex, using all these 6 C iv lines 
to calculate the associated Hi fraction leads to a biased result. 
Therefore, we used 4 C iv lines which are part of a C iv com- 
plex with less than 3 components in order to estimate the missed 
enriched H i fraction. 

There is a total of 11 Hi lines at logA^Hi = [12.75, 16.00] 
within 100 km s ' centered at the 4 weak Civ lines. The ra- 
tio of the Civ enriched Hi lines and the total Hi lines in the 
wavelength regions corresponding to the S/N > 120 Civ forest 
is [2/231, 4/315, 3/174, 1/59, 1/28, 0/10] for logA^Hi = [12.75- 
13.00, 13.0-13.5, 13.5-14.0, 14.0-14.5, 14.5-15.0, 15.0-16.0], 
respectively, or < 3%. Again the fraction of missed C iv is negli- 
gible even at the low redshift bin. 

Note that our estimate on the true undetected C iv frac- 
tion is uncertain. However, from Fig. |T3] the incompleteness at 
log A^civ = 12.2 is less than 10% or within the Icr Poisson error 

While it is clear that the incompleteness does not play a sig- 
nificant role in the H i detection down to log A^hi - 12.75 and the 
Civ detection down to logA^civ - 12.2, the combination of the 
H I and C iv detection could introduce a bias in the C iv assigning 
method. The pixel optical depth method which correlates the op- 
tical depth of H I (thi) and C iv (tciv) at the same redshift shows 
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Fig. 15. The A^Hi-A^civ diagram for log A^hi = [12.0, 17.8] from 
the high-order fit sample at 1.9 < z < 3.2. The upper panel is 
for the Avciv = ±100 km s ' sample, while the lower panel for 
the Avciv = +10 km s ' sample. Open circles represent Hi ab- 
sorbers associated with all the possible C iv components, since 
a single Hi line could be assigned to several Civ lines. On the 
other hand, red filled circles indicate a H i absorber associated 
with only one closest C iv. The total number of open (filled) cir- 
cles is 1082 (451) and 184 (163) for the Avciv = ±100 km s"' 
and Ave IV = ±10 km s"' sample, respectively. 



that at z ~ 3 there is a one-to-one positive correlation between 
the median thi and the median tciv down to logTni ~ 0.15 or 
log A^Hi ~ 13.73 for bai = 28 km s"' (a median bm of the forest 
at z ~ 2.5) dSchave et alj2003h . Below log th i ~ 0. 15, the th , sig- 
nal is blended with noise at logrciv ~ 0.001 or logA^civ < 110 
for bciv = 9.5 km s ' (a median bciv of all the C iv lines in our 
UVES sample). 

This result suggests that many Iow-A^hi absorbers might be 
mis-assigned as unenriched H i absorber in our C iv assigning 
method. Unfortunately, the lower log A^c iv ~ 1 1 limit that a typ- 
ical optical depth analysis explores is an order of magnitude 
lower than our adopted low A^civ limit of logA^civ = 12.2. This 
logA^civ ~ 11 limit cannot be obtained even in the highest S/N 
C IV region with S/N > 220 in our UVES spectra. Therefore, our 
C IV analysis can not confirm, nor refute the results from the op- 
tical depth method. 

Fig. [15] shows the Nhi-Nqiv diagram for the Avciv - 
±100kms"' sample (the upper panel) and for the Avciv - 
±10 km s ' sample (the lower panel). Since one Hi line can be 
associated with several Civ lines, data points at the same A^hi 
represent the same H i absorber Open circles show all the H i 
absorbers associated with all the possible C iv lines. Red filled 
circles indicate Hi absorbers associated with only one closest 



C IV within the search velocity range. With a larger search veloc- 
ity range, the Avciv = ±100 km s"' sample has more lines. The 
number of the red filled circles increases abruptly at log A^hi ^ 
15.0 at both redshift ranges, more prominently at the high red- 
shift bin. This is simply due to the fact that the number of weaker 
H I absorbers is larger than stronger H i absorbers. 

If our C IV assigning method were biased due to our failure 
to detect Civ lines toward lower A^hi values, there should be a 
correlation in A^hi and A^civ, such that a lower A^hi line tends to 
be associated with a lower A^c iv line (cf. the relation between 
the median thi and the median tciv) or the number of the C iv- 
enriched Hi lines at lower A^hi is smaller No such correlations 
are seen in Fig. [15] Note that our method deals with the fitted 
individual lines, while the optical depth analysis works with sta- 
tistical, median values. The optical depth analysis is not sensitive 
to any minor C iv po pulation, such as high-metallicity absorbers 
(ISchave et al.ll2007h . 

In reality, the detection of weak C iv is dependent on the lo- 
cal S/N as well as the combination of bciw and A^civ. The S/N 
of a spectrum does not change in a way to satisfy a higher S/N 
at strong Hi absorbers and a lower S/N at weaker Hi absorbers 
or vice versa. Usually the S/N changes over a larger wavelength 
interval than the wavelength interval between typical strong H i 
lines. In addition, strong and weak Hi lines do not occupy a 
portion of a spectrum separately, but exist mixed along the spec- 
trum. If a weak Civ were detected associated with a high-A^Hi 
line, a similar strength of C iv, if exists, should be detected for 
low-A^Hi lines nearby or in a similar S/N region. Therefore, un- 
less a majority C iv fraction at lower A^c iv and/or lower A^h i has a 
very large bcw value, i.e. high gas temperature, our C iv assign- 
ing method does not introduce a serious selection bias within the 
adopted A^c iv limit. 

5.2. Results 

5.2.1 . Number density evolution of tine C iv-enriclied 
absorbers 

In a similar way to the analysis of all the H i absorbers, we cal- 
culate the absorber number density evolution dnei+civ/dz on 
a quasar by quasar analysis for all the C iv-enriched H i ab- 
sorbers. The resulting d«Hi+civ/dz evolution is shown in Fig. 
[T6]forthe Avciv = ±100 km s"' and Ave iv = ±10kms"' in- 
terval from the high-order Lyman fit samples. For the Aveiv = 
±100 km s"' sample, the Ql 101-264 sightline does not show 
any C iv in the redshift range of interest due to its short red- 
shift coverage. For the Aveiv = ±10 km s"' sample, 7 sightlines 
(HE2347-4342, Q0002-422, PKS0329-255, HE1347-2457, 
Q0109-3518, Q0122-380 and Ql 101-264) out of 16 have no 
Civ-enriched Hi absorbers at logA^ei = [12.75,14.0], while 
only only sightline (Ql 101-264) has no Civ-enriched Hi ab- 
sorber at logA^Hi - [14, 17]. This is caused by the combina- 
tion of two facts that Civ tends to be associated with strong 
H I absorbers and that the small search velocity is not adequate 
due to the velocity difference between H i and C iv observed in 
many enriched absorbers. For these sightlines, dwHi+civ/dz is 0. 
Therefore, their log WHi+civ/dz is set to be with a downward 
arrow in Fig. [16] 

As for the entire Lyor forest analysis, the dnei+civ/dz evolu- 
tion resembles a power law. Therefore, linear regressions have 
been obtained from the data set and its results are summarised in 
Table [4] Sightlines showing no C iv-enriched H i absorbers were 
not included in the regression. Similar to the entire high-order- 
fit Hi sample, the Aveiv = ±100 km s"' sample shows a de- 
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Fig. 16. Quasar by quasar line number density evolution of C iv-enriched H i absorbers. The left panels are derived from the Ave iv = 
+ 100 km s ' sample, while the right panels represent the Avciv = ±10 km s ' one. The open circles in the upper panel represent 
Civ-enriched absorbers having logA^ei - [12.75, 14.0] and the filled circles in the lower panel represent log A^hi - [14, 17]. The 
vertical error bars mark 1 cr Poisson errors, while the x-axis error bars show the redshif t range covered by each sightline. Sightlines 
having no Civ-enriched Hi absorbers for a given velocity range are plotted at logddnciv/ddz = 0. The solid lines represent linear 
regressions to the data, using the parameters summarised in Table|4] Sightlines with no Civ-enriched absorbers, logdnni+civ/dz is 
plotted to be with a downward arrow. 



cline in the C iv-enriched absorber number density with decreas- 
ing redshift. This behaviour is present in both column density 
ranges. Comparing these results with the quasar-by-quasar dn/dz 
of the entire high-order fit sample at 1.9 < z < 3.2 shows that 
the Civ-enriched absorbers at logA^Hi - [14, 17] has a steeper 
slope (5.27 + 0.99), but completely consistent within Icr. The 
robust result on the d«Hi+civ/dz evolution requires a large red- 
shift coverage and more sightlines per redshift coverage, espe- 
cially at high column density range. With a lack of more C iv 
forest data at Zfoiest > 3, dnei+civ/dz derived in this study at the 
high column density range should be considered less robust com- 
pared to the entire Hi dn/dz. Similarly, the dnei+civ/dz slope 
(0.78 + 0.92) at logA^Hi = [12.75, 14.0] is also consistent with 
the one (1.38 + 0.22) of the entire high-order-fit forest sample, 
given the rather large uncertainty. The actual number densities 
are lower at both column density ranges. 

This becomes apparent in the left panel of Fig. [17] where 
the ratios of the number densities of the C iv-enriched systems 
dnHi+civ/dz and the number density of the entire sample dn/dz 
are shown. The results for the Avciv - +100 km s"' sample 
(filled circles) show that there is no significant evolution of 
the Crv enrichment fraction for logA^Hi - [12.75, 14.0]. For 
logA^Hi = [14, 17], the enrichment fraction is consistent with 
no redshift evolution, considering a large scatter at z ~ 2 and 
a lack of data at z > 3. For the low column density logA^ei = 
[12.75, 14.0] sample we find that around 5% of all the Hi ab- 
sorbers show C IV enrichment. The C iv enrichment fraction is 
higher for larger column densities of log A^hi - [14, 17], where 
around 40% of the absorbers are C iv-enriched. 

This picture changes slightly for the Avciv - ± lOkms"' 
sample. For the high column densities, the d«Hi+civ/dz evolution 



Table 4. Linear regression results for the number density evolu- 
tion dnHi+civ/dz of the Civ-enriched Hi forest absorbers in the 
quasar by quasar analysis. 





Avcv = ± 


100 km s^' 


Avciv = + 


10 km 


-1 


A log Ah, 


log A 


7 


log A 


7 




12.75 - 14.0 


0.92 ± 0.49 


0.78 ± 0.92 


2.06+1.02 


-2.49 ± 


1.94 


14.0- 17.0 


-1.40 + 0.53 


5.27 ± 0.99 


-1.22 + 0.64 


4.52 + 


1.22 



is less strong compared to the one of the Avciv - +100 km s" 
sample. However, both are still consistent within Icr due to a 
large uncertainty. Only the number density itself decreases by 
a factor of 1.7. The enrichment fractions in the right panel of 
Fig. [TTl show that now around 20% to 30% of the high column 
density H i absorbers are C iv-enriched. 

On the other hand, dnni+civ/dz increases with decreasing 
redshift for the low column densities. Its negative slope of 
y - -2.49 ± 1.94 shows an opposite behaviour from the one 
(r = 0.78 + 0.92) of the Avciv = +100 km s"' sample. This 
negative slope is in part caused by the inadequacy in our C iv 
assigning method at the small search velocity, and in part by the 
fact that the number of high-metallicity absorbers increases at 
low redshift (Schave et al. 2007). However, due to several sight- 
lines containing no C iv-enriched weak H i absorbers which are 
not included in the power-law fit, the negative slope should not 
be taken literally. The fraction of enriched absorbers increases 
from ~0.5% at z ~ 3 to ~1.5% at z ~ 2.1, as expected from 
dnHi+civ/dz at the low Hi column density. However, keep in 
mind that the cosmic variance is large as some sightlines show 
no enriched weak H i absorbers. 
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Fig. 17. The fraction of the C iv-enriched H i absorber number density to the total absorber number density as a function of redshift. 
The left panels are derived from the Avciv = +100 km s ' sample, while the right panels represent the Avciv = +10kms ' 
sample. The black open circles represent a column density interval of logA^ei = [12.75, 14.0] and the filled circles represent 
logA^Hi = [14, 17]. The vertical error bars mark Icr Poisson errors, while the x-axis error bars show the redshift range covered 
by each sightline. In the left panel for logA^Hi - [12.75, 14.0] (gray open circles), two lowest data points at log(l + z) ~ 0.5 (or 
z ~ 2) including an upper limit are from Ql 101-264 and Q0122-380. Both have a short redshift coverage, therefore become more 
susceptible to cosmic variance. Sightlines with no C iv-enriched H i absorbers are plotted as upper limits with an arbitrary value of 
0.05 and 0.002 for log A^hi = [14, 17] and log A^hi = [12.75, 14.0], respectively. 



There are two distinct groups of C iv absorbers assigned to 
the low H I column density. One group is associated with strong, 
saturated high column density H i absorbers. These absorbers are 
sometimes accompanied by lower A^hi absorbers within a veloc- 
ity range of Av < 200 km s In these systems, the C iv absorp- 
tion is usually found within 20 km s"' to the strongest Hi lines 
(Kim et al. 2013, in preparation). Therefore, these accompanied 
low H I column density systems get associated with the C iv ab- 
sorbers if the velocity range Ave iv is large. With a small veloc- 
ity search interval, however, only H i systems that have C iv in 
their direct vicinity are flagged as C iv-enriched. This means that 
the aforementioned low column density systems around strong 
absorbers are not considered C iv-enriched in a small velocity 
search interval. 

Another C iv-enriched group consists of usually isolated, low 
column density H i absorbers associated with strong C iv ab- 
sor ption, i.e. the same high-metallicity forest population studied 
by ISchave et alj (l2007i) . An example of such a system toward 
HEl 122-1648 is shown in Fig. [18] In this velocity plot, an Hi 
absorption feature is hardly recognisable, while strong C iv and 
N V doublets are present. The existence of both doublets makes 
the identification of this absorber secure. Due to the low A^hi 
and high A^metais, these systems show a higher ionisation and 
a higher metallicity compared to a typical absor b er with simi- 
lar A^h i (ICarsweU et al."2002t [Schave et al.ll200"7l) . ISchave et al.l 
(|2QQ7|) speculate that these systems could be responsible for 
transporting metals from galaxies to the surrounding IGM. As 
the velocity difference between Hi and metal lines for these 
systems are usually very small, they dominate the weaker C iv- 
enriched forest at log A^hi < 14 for the Avciv - ±10 km s"'. In 
addition, the high-metallicity absorbers are more common at low 
redshift. 

The different characteristics of these two C iv groups ex- 
plains the different dnciv/dz behaviour between the Avciv - 
+ 100 kms"' and Avciv = ±10kms"' samples at logA^Hi - 
[12.75, 14.0]. With recent observational evidence that metals are 
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Fig. 18. A velocity plot of a highly enriched C iv absorber at z = 
2.030069 toward HEl 122-1648. The zero velocity is centered at 
z = 2.030069. Although there is no obvious Lya absorption seen 
at the zero velocity, both C iv and N v doublets are present to 
secure the existence of this absorber Note that the y-axis range 
for each ion is different: from the normalised flux to 1 for H i 
and from 0.5 to 1 for the rest of the ions. 
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Fig. 19. The distribution function for C iv-enriched H i lines for 
Avciv = ±100kms~^ (red filled squares) and for Avciv - 
+ 10kms"' (blue stars) at 1.9 < z < 3.2. Also shown is 
the differential column density distribution function for all H i 
Lya absorbers excluding Q0055-269 and J2233-606 (black 
filled circles) in the same redshift range analysed for the C iv- 
enriched forest. The solid line indicates the fit to filled circles 
forlogA^Hi = [12.75, 14.0] : logdN/idNmdX) = (7.43+0.44) + 
(-1.44 + 0. 03)xlogA^Hi- The vertical errors indicate IcrPoisson 
errors, while the x-axis error bars show the A^hi range covered 
by each data point. All the grey data points indicate that the data 
are incomplete at log A^hi < 12.75. 



only found close to galaxies in the circ um-galactic medium at 
2 < z < 4 and not far away from galaxies (" Adelberger et al.lF2005t 
ISteidel et al.i r2010), our results could provide further theoreti- 
cal constraints for this interpretation. It could well be that the 
high-metallicity forest population is completely different from 
the typical, low-metallicity forest and resides in a different inter- 
galactic space. However, due to the low gas density and high 
temperature, the Lya forest does not have in-situ star forma- 
tion. Metals associated with the Hi forest should have been 
transported from nearby galaxies. In other words, all the C iv- 
enriched absorbers are close to galaxies. 

5.2.2. Differential column density distribution function of the 
C iv-enriched forest 

Fig. [19] shows the differential column density distribution func- 
tion for C iv-enriched H i absorbers for 1 .9 < z < 3.2. Red filled 
squares and blue stars represent the search velocity ranges of 



Table 5. Linear regression results for the differential column 
density distribution of the Civ-enriched forest at logA^Hi = 
[14.5, 17.0]. Here the normalisation point log (d«/(dA^Hi d^))o is 
denoted by B. 
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0.21 
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6.90+1.45 


-1.44 + 0.10 


5.97 + 1.94 
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0.13 



Avmetai = ±100 km s^' and Avmetai = +10 km s"', respectively. 
Black filled circles are for all H i lines (excluding J2233-606 and 
Q0055-269), regardless of their metal association. As in Figs|7] 
and[TOl the binsize of log A^H I = 0.25 is used at log A^H I = [ 1 2 , 1 5 ] , 
then the binsize of 0.5 at logA^Hi = [15,18]. The solid line 
indicates the fit to filled circles for logA^Hi = [12.75, 14.0] : 
\ogdNliANnAX) = (7.43 + 0.44) + (-1.44 + 0.03) x logA^Hi- 
The total absorption distance is Xiz) - 19.6652 for the redshift 
ranges analysed in this subsection. 

For logA^Hi > 15, the CDDF of the enriched forest is not 
sensitive to our choice of the search velocity and the Avmetai = 
+ 100 km s ' CDDF becomes almost identical with the CDDF 
of the entire Hi sample. For the column densities logA^Hi = 
[14, 17], the CDDF functional form of the enriched forest shows 
a power-law with a similar slope obtained for the entire H i ab- 
sorbers at logA^Hi - [12.75, 14.0], but with a smaller normalisa- 
tion value. 

At log A^H I < 15, the distribution function of the C iv-enriched 
forest starts to deviate significantly from the CDDF of the 
entire Hi sample. The CDDF of the Civ-enriched Hi forest 
starts to flatten out toward lower A^hi at both search velocity 
ranges. Furthermore the flattening of the enriched forest depends 
strongly on the choice of Avciv The large search velocity results 
in a steeper slope with a less fluctuation than the small one. This 
is due to the Avciv - ±10 km s"' sample being predominantly 
sensitive to highly enriched absorbers at logA^Hi < 14 and less 
sensitive to mis-aligned broad Civ lines with b > \Q km s"'. 
Note that our method to associate H i with C iv is only depen- 
dent on the relative velocity difference between the line centers, 
but not the C iv profile shape. The large velocity range includes 
broader C iv lines up to ~ 100 km s ' as well as narrow, highly 
enriched absorbers. The Avciv - ±100 km s"' velocity range is 
a better filter to associate H i and C iv. 

The flattening of the distribution function seen at logA^Hi < 
15 by Civ-enriched absorbers cannot be caused by the incom- 
pleteness of the H I sample. The H i incompleteness would result 
in a similar flattening as is seen at logA^Hi < 12.75 for the entire 
sample (as seen in Fig.|7]i. However, our sample of H i absorbers 
is complete for column densities larger than log A^h i > 12.75. 

As discussed in Section 15.11 the flattening of the Civ- 
enriched forest at logA^Hi < 15 could be in part caused by the 
missed weak, broad C iv lines. However, Fig. [13] shows that the 
C IV CDDF at log A^c iv > 12.2 is not strongly affected by the C iv 
incompleteness. The number ratio of the entire H i forest lines 
and the Civ-enriched forest lines at logA^Hi ~ 13 is ~ 25 for 
the Avciv - ±100 km s"' sample. This ratio increases to ~ 260 
for the Ave IV = ±10 km s"' sample. Even if we took a max- 
imum correction for the C iv incompleteness of 50%, roughly 
consistent with the results by Giallongo et al. ( 1996) (see their 
Section 2.3), the CDDF flattening of the enriched forest toward 
lower A^Hi is still present. Therefore, this flattening is real and 
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Fig. 20. Differential column density distribution function for enriched and unenriched absorbers in our high-order fit sample us- 
ing Avciv - ±100 km s"' (upper panels) and Avciv = +10kms"' (lower panels) at the two different redshift ranges. Black 
filled circles and gray stars mark unenriched absorbers. Red filled squares and blue stars indicate the C iv-enriched forest for the 
Ave IV - +100 km s"' sample and the Avciv = +10 km s"' sample, respectively. Both Q0055-269 and J2233-606 are excluded in 
the analysis. The solid black line indicates the fit to the entire Hi sample and the whole redshift range at log A^hi = [12.75, 14.0] as 
in Fig. [19] Red and blue dashed lines represent the fit to each C iv-enriched forest sample at log A^hi = [14.5, 17]. The vertical errors 
indicate Icr Poisson errors, while the x-axis error bars show the A^hi range covered by each data point. The lower parts of the panels 
show the difference between the observed CDDF and the expected CD DP from the power-law fit obtained for the entire H i samples 
(black sohd hnes). 



physically related to the characteristics of the C iv-enriched ab- 
sorbers only with A^civ > 12.2. 

The observation that the differential column density distribu- 
tion of the C iv-enriched forest flattens at low column densities 
can be easily explained by the fact that the enrichment fraction 
with log A^civ > 12.2 becomes smaller as A^hi decreases. 



At logA^Hi = [13.0, 13.5, 14.0, 14.5, 15.0, 15.5], the fraction 
of the metal enriched forest for the Ave iv = ± 1 00 km s" ' sample 
is roughly [4,6, 11,31,49,78]%, respectively. This enrichment 
fraction can be roughly inferred from the difference between the 
entire H i CDDF and C iv-enriched CDDF in Fig. [19] 
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The different CDDF shape between the C iv-enriched ab- 
sorbers and unenriched absorbers strongly supports that the C iv- 
enriched absorbers arise from the different physical environ- 
ment, i.e. the circum-galactic medium, while the unenriched for- 
est has its origin as the intergalactic medium. The fact that the 
number of Civ-enriched absorbers decreases with decreasing 
A^Hi is also consistent with the picture of IGM metal enrich- 
ment models by galactic winds (Aguirre et al. 2001). The lower 
the H I column density of absorbers is, the farther they are from 
high-density gas concentrations where galaxies are formed. As 
galactic winds have a limited life time and outflow velocity to 
transport metals in to the low-density IGM, weaker absorbers 
will not be likely to be metal enriched. 

The redshift evolution of the distribution function of C iv- 
enriched and unenriched absorbers (not the entire H i absorbers) 
is shown in the upper panels of Fig. |20] for the Avciv - 
+ 100 km s ' sample and in the lower panels for the Avciv - 
ilOkms"' sample at the redshift ranges z = [1.9,2.4] and 
[2.4, 3.2]. The total absorption distance is X(z) = 11.0940 and 
8.57061 at z = [1.9,2.4] and [2.4, 3.2], respectively, exclud- 
ing Q0055-269 and J2233-606. To increase the absorber num- 
ber for each A^hi bin, the binsize of logA^Hi - 0.4 is used at 
log A^Hi = [12.75, 17.95]. Black filled circles and gray stars mark 
absorbers without C iv. Red filled squares and blue stars indi- 
cate the C iv-enriched forest for the Av'civ = ±100 km s"' sam- 
ple and the Avciv = ±10 km s"' sample, respectively. The solid 
black line indicates the power-law fit to the entire H i sample at 
1.9 < z < 3.2 and at logA^Hi = [12.75, 14.0] as in Fig. \TU\ Red 
and blue dashed lines represent the power-law fit to each C iv- 
enriched forest sample at log A^hi - [14, 17] (see Table|5]). 

The upper panel of Fig. |20]suggests that the entire absorber 
population can be considered as the combination of two pop- 
ulations of well-characterised absorbers, the enriched absorbers 
and the unenriched absorbers. The C iv-enriched absorbers dom- 
inate at log A^Hi > 15. Their CDDF is well-described as a power 

law. The slope /3 1.45 obtained at logA^Hi = [14.5, 17] (red 

dashed lines) is similar to the slope /3 ~ -1 .44 for the entire H i 
sample atlogA^ei - [12.75,14.0] at both redshifts. The normali- 
sation value for the C iv-enriched forest is smaller, with about 10 
times lower absorber numbers. The enriched absorbers do not 
show any strong redshift evolution at logA^Hi > 15, while the 
CDDF flattening at logA^Hi < 15 seems to be weaker at the low 
redshift. The unenriched absorbers dominate at log A^hi < 15 with 
a power-law CDDF. At higher A^hi, the unenriched absorbers 
become significantly deviated from the extrapolated power law 
obtained at lower A^hi- There are no unenriched absorbers at 
logA^Hi>16.0. 

The lower panels for the Avciv = ±10 km s ' sample show 
similar results. There is no strong redshift evolution for the C iv- 
enriched absorbers at logA^Hi > 15. Again, there is a sugges- 
tion that the flattening at logA^Hi < 15 becomes less significant 
at lower redshifts. The unenriched forest starts to dominate at 
log A^Hi < 15. However, the Avciv = ±10 km s"' sample shows a 
less-smooth CDDF. The highest-A^ei data point at z = [2.4, 3.2] 
illustrates the inadequacy of our assigning method of C iv to H i 
when a small search velocity was used. Absorbers contributing 
this data point are part of multi-component high column density 
systems. Their associated C iv shows a rather simple, broad, but 
mis-aligned profile from the H i center of saturated Lya profiles. 
The H I line center of some components resolved at high-order 
Lyman lines is sometimes at > 10 km s"' from the closest Civ 
component and thus they are flagged as the unenriched forest. 

In both samples, the two populations overlap at log A^hi ~ 15. 
The deviation from the power law starting at log A^hi ~ 14.0-14.5 



Table 6. Hi density relative to the critical density at logA^Hi = 
[12.75,17.0]: the entire Hi forest Qhi, the enriched forest 
^^Hi+civ and the unenriched forest Qhi-civ- Units are lO^^/i"'. 
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shown in the CDDF for the entire H i sample in Fig.[TO]is a result 
of combining two different populations which show a different 
CDDF shape. 

To obtain a rough idea on the Hi density relative to the 
critical density of the entire forest (Qhi), the enriched for- 
est (Qhi+civ) and the unen riched forest (Qhi-Civ), we used 
Eq. (A12) from'Schaye'(^l) at logA^Hi = [12.75, 17.0], by di- 
rectly integrating the observed CDDF. The ionisation fraction of 
H I to obtain the total hydrogen density Oh was not corrected, 
since it is highly uncertain. The resulting mass fractions for 
log A^Hi = [12.75, 17.0] are given in Tabled The derived CDDF 
is not well-constrained at logA^Hi ^16 due to the low number 
statistics. Therefore, the derived O values in Table |6] are only 
rough numbers. The ratio of Ohi+civ and Qhi might decrease by 
a factor of 2 at low redshift for both Avciv samples. However, a 
high uncertainty in deriving Q values does not allow to assure 
this decrease. The C iv-enriched forest accounts for ~ 40% of 
the entire forest in mass at 1.9 < z < 3.2 for the Avciv = 100 
sample. 

6. Conclusions 

Based on an in-depth Voigt profile fitting analysis of 18 high- 
redshift quasars obtained from the ESO VLT/UVES archive, we 
have studied ~ 3100 Hi absorbers to investigate the number 
density evolution and the differential column density distribu- 
tion function at 1.9 < z < 3.2 and for logA^ei = [12.75, 17.0]. 
Two methods of the Voigt profile fitting analysis have been ap- 
plied, one by fitting absorption profiles only to the Lyo' transition 
and another by including higher order Lyman transitions such as 
Ly/3 and Lyy. These higher order transitions provide a more re- 
liable column density measurement of saturated absorption sys- 
tems, since saturated and blended lines often become unsatu- 
rated at higher order transitions. This also enables us to resolve 
the structure of absorbers more reliably. This study has increased 
the sample size by a factor of 3 from previous similar studies at 
z > 2. In addition, we have investigated whether there exist any 
differences in the A^hi evolution between the C iv-enriched forest 
and the unenriched forest. 

We have found that the results based on the Lya-only fit are 
in good agreement with previous results on a quasar by quasar 
analysis. For our data only (values in parenthesis indicate re- 
sults including high-quality data from the literature at z > 1), the 
number density dn/dzis dn/dz = (1.46 ± 0.1 1) x (1 +2)i-67±o.2i 
((1.52 ± 0.05) X (1 -H 2)i-5i±('(») ^ (0.03 ± 0.20) x 

(1 + 2)3-40±0.36 ((Q 72 + 0.08) X (1 H- 2)216±014) at log A^hi = [13.1, 

14.0] and [14, 17], respectively. The noticeable difference in the 
exponent between our sample and the sample including the data 
from the literature for stronger absorbers is caused by the fact 
that our sample does not cover a large redshift range and that 
the evolution of dn/dz is more significant for stronger absorbers. 
The scatter between different sightlines becomes larger at lower 



21 



T.-S. Kim et al.: The high-redshift Hi and C iv forest 



redshifts and stronger absorbers due to the evolution of the large- 
scale structure. 

Combining our Lya-only fit analysis at 1.9 < z < 3.6 with 
the high-quality literature data at 0.0 < z < 4, the mean number 
density evolution is not well described by a single power law 
and strongly suggests that its evolution slows down at z < 1 .5 
at both high and low column density ranges. Although a single 
power law does not give a good description, the number density 

is d«/dz oc (1 +2)0.89±0.06 ^jjjj ^ +^y.61±0.12 log A?H, = 

[13.1, 14.0] and [14, 17], respectively. 

The differential column density distribution function 
(CDDF) from the Lya-only fit analysis is also consistent with 
previous results. The single power-law exponent is -1.44 + 0.02 
at 1 .9 <z < 3.2 and log A^hi = [12.75-14.0], with a deviation from 
the power law at log A^hi > 14.0 - 14.5. 

The high-order Lyman fits do not show any significantly dif- 
ferent results from the ones based on the Lya-only fits. The 
dn/dz evolution based on a quasar by quasar analysis yields a 
very similar result to the Lya-only fit. The mean dn/dz based 
on the combined sample from our quasars at 1 .9 < z < 3.2 is 
dn/dz = (1.89 + 0.13) (1 + 2)i-28±0-24 and dn/dz = (-0.65 + 
0.36) (1 + 2)4«±o.66 logA^jjj = [12.75, 14.0] and [14, 17], re- 
spectively. 

Using the high-order fits, we have derived the differential 
column density distribution function at 1.9 <z<3.2 and confirm 
the existence of a dip at logA^ei = [14, 18] as seen in the Lya- 
only-fit CDDF analysis. At 1 .9 <z < 3.2, the power-law exponent 
of the differential column density distribution function is - 1 .44 + 
0.03, -1.67 + 0.09 and -1.55 + 0.08 at logA^Hi = [12.75, 14.0], 
[14, 15] and [15, 18], respectively. 

By obtaining the differential column density distribution 
function for two redshift bins z = [1.9,2.4], and [2.4, 3.2], 
we observe that a deviation from the expected power-law at 
logA^Hi - [14.0, 18.0] is more prominent at lower redshifts. In 
addition, the power-law seems to be slightly steeper at the low 
redshift for the column density range log A^hi = [12.75, 14.0] in 
which the distribution function follows a perfect single power 
law. However, the CDDF at two redshift bins is consistent with 
no redshift evolution within 2cr. 

Further, we have split the entire Hi absorbers excluding 2 
quasars with a lower S/N Civ region into two samples: ab- 
sorbers associated with Civ tracing the metal enriched forest, 
and absorbers associated with no C iv tracing the unenriched 
forest. A H i absorber is considered C iv-enriched, if a Civ line 
with logA^civ greater than a threshold value is found within a 
given search velocity interval centered at each Hi absorption 
center. The threshold logA^civ = 12.2 was used since the Civ 
distribution function and the Ncw-bciv diagram show that the 
Civ detection is reasonably complete down to logA^civ - 12.2 
for a typical bcw value found at logA^civ ^ 12.2 in our sam- 
ple. We used two arbitrarily chosen search velocity intervals, 
Ave IV = +100kms"' and Avciv = +10kms"'. 

At logA^Hi = [14, 17], the dnni+civ/dz of the Civ-enriched 
Hi absorbers show a similar evolution compared to the one 
of the entire Lya forest, with a power-law decrease in num- 
ber density with decreasing redshift. The power-law slope is 
[0.78+0.92,5.27 + 0.99] forlogA^ei = [12.75,14.0] and [14, 
17] at 1.9 < z < 3.2 for the Avciv = ± 100 km s"' sample. 

The enriched fraction is fairly constant with redshift at 1 .9 < 
z<3.2. About 5% of all absorbers show an association with C iv 
at log A^Hi = [12.75, 14], while about 40% are metal enriched at 
log A^Hi = [14, 17] for the Avciv = ± 100 km s"' sample. 

For Ave IV - ± 10 km s"' sample, the low column density 
enriched absorber suggests that dnni+civ/dz increases as red- 



shift decreases, i.e. a negative slope of -2.49 + 1.94. Part of this 
behaviour is caused by the fact that high-metallicity absorbers 
which are more sensitive to the small search velocity become 
more abundant at low redshift. However, this negative evolution 
should not be taken literally since about a half of sightlines does 
not show enriched absorbers at log A^hi = [12.75, 14.0]. 

The differential column density distribution function for the 
enriched and unenriched systems show a significant difference. 
However, each shows a well-characterised CDDF. At log A^hi ^ 
15.0, the unenriched forest dominates and its distribution shows 
a power law similar to the entire forest sample. On the other 
hand, the C iv-enriched forest is found to flatten out at log A^hi ^ 
15. Depending on the search velocity interval, the number of en- 
riched systems is a factor of 25 (Avciv = +100 km s"') to 260 
(Ave IV = ±10 km s"') lower than the one of the unenriched sys- 
tems at log A^Hi = 13. This flattening is mainly caused by the fact 
that the enriched fraction of the Lya forest decreases as log A^h i 
decreases. 

At the higher A^hi range, the Civ-enriched forest domi- 
nates. Its distribution function can be described as a power law 
with its slope of -1.45 + 0.08 similar to the power-law slope 
(-1.44 + 0.03) of the entire H i forest at log A^hi = [12.75, 14.0], 
but a lower normalisation value, i.e. ~ 10 times lower in the ab- 
sorber number. The unenriched forest disappears very rapidly as 
log A^Hi increases. 

The distribution function of the entire H i forest can be de- 
scribed as the combination of these two well-characterised pop- 
ulations, overlapping at logA^Hi ~ 15. The deviation from the 
power law at logA^Hi = [14, 17] seen in the CDDF for the en- 
tire H I sample is a result of combining two different H i popu- 
lations with a different CDDF shape. This result supports other 
observational evidence from absorber-galaxy studies at z ~ 3, 
namely that metals associated with the high-redshift Lya for- 
est are within ^100 kpc of galaxies (Adelbergeret al. 20051 
Steidel e t al.l l2010i) . Absorber-galaxy studies suggest that the 
C iv-enriched and unenriched forest would arise from the dif- 
ferent spatial and physical locations, therefore having a differ- 
ent physical/evolutionary behaviour suggested by the different 
CDDF shape. Therefore, our results combined with absorber- 
galaxy studies indicate that the Civ-enriched forest is the 
circum-galactic medium, while the unenriched forest has its ori- 
gin as the intergalactic medium. 

At 1.9 < z < 3.2, the Civ-enriched forest contributes ~40% 
of the entire forest mass to the H i density relative to the critical 
density for the Aveiv = ± 100 km s"' sample. 
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Appendix A: Data 

In this appendix we present all the data for the quasar by quasar 
number density evolution (Table IATTT i. the mean number density 
evolution (Table IA.2b . the differential column density distribu- 
tion (CDDF) of the entire H i sample (Table lA3Tl, the CDDF of 
the C iv-enriched forest (Table I A. 4) and of the unenriched for- 
est (Table IA3Ti for Avwai = +100 km s"', and the CDDF of 
the C iv-enriched forest (Table lAT6t and of the unenriched forest 
(Table lAJli for Avmetai = ±10 km s"' . 
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Table A.l. Number density evolution data for each quasar 



Lya-only fit high-order fit 

logA^H. = [12.75, 14.0] logA^Hi = [14, 17] logA^H, = [12.75, 14.0] log^H, = [14, 17] 







A/ 






< ^ > 


A/ 
^z 






Q0055-269 


3.270 


0.669 


2.672±0.024 


2.157±0.042 


3.071 


0.269 


2.600±0.040 


2.269±0.057 


PKS2126-158 


3.010 


0.390 


2.642±0.032 


1.886±0.073 


3.010 


0.390 


2.642±0.032 


1.965+0.067 


Q0420-388 


2.759 


0.558 


2.656+0.026 


2.117±0.048 


2.851 


0.373 


2.722 + 0.022 


2.160+0.055 


HE0940-1050 


2.729 


0.554 


2.628+0.028 


1.997±0.055 


2.729 


0.554 


2.635±0.027 


2.070+0.051 


HE2347-4342 


2.577 


0.484 


2.602±0.030 


1.918+0.064 


2.577 


0.484 


2.600+0.030 


1.918±0.064 


Q0002-422 


2.457 


0.496 


2.651+0.029 


1.948±0.066 


2.457 


0.496 


2.659+0.028 


1.995+0.058 


PKS0329-255 


2.395 


0.513 


2.483+0.035 


1.846±0.073 


2.395 


0.513 


2.499±0.033 


1.781+0.072 


00453-423" 


2.340 


0.354 


2.665±0.033 


1.812±0.082 


2.340 


0.354 


2.663±0.033 


1.913±0.074 


HE 1347-2457 


2.300 


0.505 


2.506+0.033 


1.658±0.082 


2.300 


0.505 


2.525±0.032 


1.695+0.079 


00329-385 


2.139 


0.475 


2.571+0.031 


1.602±0.090 


2.139 


0.475 


2.569±0.032 


1.666+0.084 


HE2217-2818 


2.126 


0.479 


2.575+0.031 


1.524±0.097 


2.167 


0.396 


2.587±0.034 


1.607+0.097 


Q0109-3518 


2.127 


0.443 


2.518+0.035 


1.716±0.072 


2.158 


0.380 


2.507±0.038 


1.721+0.088 


HE 1122-1648 


2.124 


0.467 


2.588+0.031 


1.793±0.074 


2.129 


0.458 


2.589±0.031 


1.816±0.073 


J2233-606 


1.978 


0.445 


2.533±0.034 


1.582+0.094 


2.086 


0.230 


2.501+0.048 


1.416±0.149 


PKS0237-23 


1.972 


0.415 


2.519+0.036 


1.613±0.094 


2.070 


0.219 


2.492±0.050 


1.505+0.139 


PKS 1448-232 


1.947 


0.456 


2.514+0.034 


1.773±0.076 


2.064 


0.222 


2.517±0.048 


1.800+0.103 


Q0122-380 


1.921 


0.442 


2.434+0.038 


1.735±0.081 


2.059 


0.165 


2.456±0.059 


1.687+0.131 


Ql 101-264 


1.989 


0.216 


2.493+0.050 


1.512±0.139 


1.998 


0.197 


2.460±0.054 


1.307+0.176 


" It includes a 


sub-DLA, which introduces a gap in the Lya redshift range. The calculation of dn/dz and the Usted Az take into account the 


redshift gap. 


















Table A.2. Mean number density evolution data for Az = 0.26 
















Lya-only fit 




high-order fit 








log 


AfHi = [12.75, 14.0] 


log A'h, = [14,17] loi 


;iVH, = [12.75, 14.0] 


logA^H, = [14, 17] 






< z >" 




log AnIAz 


logd;j/dz 


log dn/dz 


log dn/dz 






2.067 




2.543±0.016 


1.630+0.044 


2.540+0.0160 


1.626±0.0442 






2.263 




2.516±0.017 


1.680±0.044 


2.542±0.0173 


1.723±0.0422 






2.552 




2.542+0.020 


1.874±0.042 


2.545±0.0199 


1.883±0.0416 






2.813 




2.647+0.022 


2.024±0.045 


2.664±0.0219 


2.094+0.0412 






3.027 




2.640±0.024 


2.078±0.045 


2.650±0.0240 


2.136±0.0424 





The exactly same redshift range was used for both fit samples. 



Table A.3. CDDF of the entire H i forest / = log (dA?/(dA?Hi dX)) for the high-order fit sample 





z = 


1.9-3.2 




z = 


1.9-2.4 




z = 


2.4-3.2 




logA^Hl 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


12.875 


-11.077 


0.015 


0.016 


-11.082 


0.021 


0.022 


-11.071 


0.022 


0.023 


13.125 


-11.437 


0.017 


0.018 


-11.441 


0.023 


0.024 


-11.432 


0.025 


0.026 


13.375 


-11.795 


0.019 


0.020 


-11.840 


0.027 


0.029 


-11.747 


0.027 


0.028 


13.625 


-12.163 


0.022 


0.023 


-12.206 


0.031 


0.034 


-12.117 


0.031 


0.033 


13.875 


-12.507 


0.024 


0.026 


-12.555 


0.035 


0.038 


-12.457 


0.034 


0.037 


14.125 


-12.992 


0.032 


0.034 


-13.118 


0.049 


0.055 


-12.879 


0.041 


0.045 


14.375 


-13.459 


0.040 


0.045 


-13.611 


0.064 


0.075 


-13.329 


0.051 


0.058 


14.625 


-13.781 


0.044 


0.049 


-14.083 


0.081 


0.099 


-13.579 


0.051 


0.058 


14.875 


-14.290 


0.058 


0.067 


-14.351 


0.082 


0.102 


-14.227 


0.078 


0.095 


15.250 


-14.933 


0.055 


0.063 


-15.203 


0.097 


0.125 


-14.744 


0.065 


0.076 


15.750 


-15.970 


0.097 


0.125 


-16.129 


0.149 


0.228 


-15.835 


0.119 


0.165 


16.250 


-16.674 


0.119 


0.165 


-16.629 


0.149 


0.228 


-16.733 


0.176 


0.301 


16.750 


-17.219 


0.125 


0.176 


-17.907 


0.301 




-16.932 


0.131 


0.189 


17.250 


-17.896 


0.149 


0.228 


-18.106 


0.232 


0.533 


-17.733 


0.176 


0.301 


17.750 


-19.174 


0.301 










-18.835 


0.301 
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Table A.4. CDDF / = log {dN/idNm dX)) of the C iv-enriched forest for Avmetai = ±100 km s"' 





z - 


1.9-3.2 




Z - 


1.9-2.4 




Z = 


2.4-3.2 




log A'h I 


f 


+A,f 


-A,f 


f 


+Ay 


-A/ 


f 


+A/ 


-A/ 


12.950 


-12.655 


0.066 


0.078 


-12.632 


0.084 


0.104 


-12.686 


0.100 


0.130 


13.350 


-13.021 


0.064 


0.075 


-13.096 


0.090 


0.113 


-12.940 


0.086 


0.107 


13.750 


-13.410 


0.063 


0.074 


-13.519 


0.092 


0.117 


-13.301 


0.082 


0.102 


14.150 


-13.855 


0.066 


0.078 


-13.944 


0.094 


0.121 


-13.762 


0.088 


0.110 


14.550 


-14.232 


0.065 


0.076 


-14.495 


0.110 


0.148 


-14.031 


0.076 


0.093 


14.950 


-14.705 


0.070 


0.083 


-14.861 


0.106 


0.141 


-14.562 


0.088 


0.110 


15.350 


-15.344 


0.090 


0.113 


-15.530 


0.139 


0.206 


-15.183 


0.110 


0.148 


15.750 


-16.069 


0.125 


0.176 


-16.173 


0.176 


0.301 


-15.964 


0.161 


0.257 


16.150 


-16.645 


0.149 


0.228 


-16.697 


0.198 


0.374 


-16.585 


0.198 


0.374 


16.550 


-16.978 


0.139 


0.206 


-17.274 


0.232 


0.533 


-16.764 


0.161 


0.257 


16.950 


-17.746 


0.198 


0.374 








-17.385 


0.198 


0.374 


17.350 


-18.146 


0.198 


0.374 


-18.074 


0.232 


0.533 


-18.263 


0.301 


oo 


17.750 


-19.023 


0.301 


oo 






oo 


-18.663 


0.301 


oo 



Table A.5. CDDF / = log {dN/idNm dX)) of the unenriched forest for Avmetai = ±100 km s"^ 



z = 1.9 - 3.2 z = 1.9 - 2.4 z = 2.4 - 3.2 



lOgAfHl 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


12.950 


-11.200 


0.013 


0.014 


-11.210 


0.018 


0.018 


-11.187 


0.020 


0.020 


13.350 


-11.797 


0.016 


0.017 


-11.831 


0.023 


0.024 


-11.757 


0.024 


0.025 


13.750 


-12.374 


0.020 


0.021 


-12.389 


0.027 


0.029 


-12.355 


0.029 


0.032 


14.150 


-13.114 


0.029 


0.032 


-13.240 


0.044 


0.050 


-12.991 


0.038 


0.042 


14.550 


-13.920 


0.046 


0.051 


-14.097 


0.073 


0.088 


-13.764 


0.057 


0.066 


14.950 


-14.718 


0.071 


0.084 


-14.828 


0.103 


0.135 


-14.607 


0.092 


0.117 


15.350 


-15.845 


0.149 


0.228 


-15.773 


0.176 


0.301 


-15.962 


0.232 


0.533 


15.750 


-16.722 


0.232 


0.533 


-16.775 


0.301 


oo 


-16.663 


0.301 


00 


16.150 


-17.122 


0.232 


0.533 


-16.874 


0.232 


0.533 









Table A.6. CDDF / = log {dN/idNui dX)) of the C rv-enriched forest for Avmetai = ±10 km s"' 



z= 1.9 -3.2 z= 1.9 -2.4 z = 2.4 -3.2 



logAfRi 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


/ 


+A/ 


-A/ 


12.950 


-13.524 


0.161 


0.257 


-13.276 


0.161 


0.257 








13.350 


-14.146 


0.198 


0.374 


-13.897 


0.198 


0.374 








13.750 


-14.023 


0.119 


0.165 


-14.076 


0.161 


0.257 


-13.964 


0.161 


0.257 


14.150 


-14.277 


0.103 


0.135 


-14.220 


0.125 


0.176 


-14.364 


0.161 


0.257 


14.550 


-14.544 


0.090 


0.113 


-14.973 


0.176 


0.301 


-14.286 


0.100 


0.130 


14.950 


-14.861 


0.082 


0.102 


-15.020 


0.125 


0.176 


-14.716 


0.103 


0.135 


15.350 


-15.477 


0.103 


0.135 


-15.773 


0.176 


0.301 


-15.263 


0.119 


0.165 


15.750 


-16.178 


0.139 


0.206 


-16.297 


0.198 


0.374 


-16.060 


0.176 


0.301 


16.150 


-16.645 


0.149 


0.228 


-16.697 


0.198 


0.374 


-16.585 


0.198 


0.374 


16.550 


-17.124 


0.161 


0.257 


-17.274 


0.232 


0.533 


-16.985 


0.198 


0.374 


16.950 


-17.922 


0.232 


0.533 








-17.561 


0.232 


0.533 


17.350 


-18.623 


0.301 


oo 








-18.263 


0.301 


00 
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Table A.7. CDDF / = log {dN/idNm dX)) of the unenriched forest for Avmetai = ±10 km s"^ 



z= 1.9 -3.2 z= 1.9 -2.4 z = 2.4 -3.2 



1 n o (\', , 


f 
J 




-A f 












f 
1 


+A f 


-A f 


12.950 


-11.187 


0.013 


0.013 


-11 


.197 


0.017 


0.018 


-11, 


,173 


0.019 


0.020 


13.350 


-nil A 


0.016 


0.017 


-11 


.811 


0.022 


0.023 


-11, 


,730 


0.023 


0.024 


13.750 


-12.345 


0.019 


0.020 


-12, 


.366 


0.026 


0.028 


-12, 


,318 


0.028 


0.030 


14.150 


-13.067 


0.028 


0.030 


-13, 


.201 


0.043 


0.047 


-12, 


,939 


0.036 


0.039 


14.550 


-13.823 


0.041 


0.046 


-13, 


.995 


0.065 


0.077 


-13, 


,670 


0.052 


0.059 


14.950 


-14.600 


0.062 


0.073 


-14, 


.719 


0.092 


0.117 


-14, 


.482 


0.081 


0.099 


15.350 


-15.582 


0.114 


0.156 


-15, 


.530 


0.139 


0.206 


-15, 


.660 


0.176 


0.301 


15.750 


-16.421 


0.176 


0.301 


-16, 


.474 


0.232 


0.533 


-16, 


.361 


0.232 


0.533 


16.150 


-17.122 


0.232 


0.533 


-16, 


.874 


0.232 


0.533 










16.550 


-17.522 


0.232 


0.533 










-17, 


,161 


0.232 


0.533 


16.950 


-18.223 


0.301 


CO 










-17, 


,863 


0.301 


OO 


17.350 


-18.322 


0.232 


0.533 


-18, 


.074 


0.232 


0.533 










17.750 


-19.023 


0.301 


oo 










-18, 


.663 


0.301 


OO 
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